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The generation of a baryon asymmetry via leptogenesis is usually studied by means of clas- 
sical kinetic equations whose applicability to processes in the hot and expanding early universe 
is questionable. The approximations implied by the state-of-the-art description can be tested in 
a first-principle approach based on nonequilibrium field theory techniques. Here, we apply the 
Schwinger-Keldysh/Kadanoff-Baym formalism to a simple toy model of leptogenesis. We find that, 
within the toy model, medium effects increase the vertex contribution to the CP-violating param- 
eter. At high temperatures it is a few times larger than in vacuum and asymptotically reaches 
the vacuum value as the temperature decreases. Contrary to the results obtained earlier in the 
framework of thermal field theory, the corrections are only linear in the particle number densities. 
An important feature of the Kadanoff-Baym formalism is that it is free of the double-counting 
problem, i.e. no need for real intermediate state subtraction arises. In particular, this means that 
the structure of the equations automatically ensures that the asymmetry vanishes in equilibrium. 
These results give a first glimpse into a number of new and interesting effects that can be studied 
in the framework of nonequilibrium field theory. 

PACS numbers: ll.10.Wx, 98.80.Cq 
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I. INTRODUCTION 

The almost complete absence of antimatter on Earth, 
in the solar system and in hadronic cosmic rays suggests 
that the universe is baryonically asymmetric. This con- 
clusion is confirmed by experimental data on the abun- 
dances of the light elements [l[ and precise measurements 
of the cosmic microwave background spectrum 0, H[ . 

The baryon asymmetry of the universe can be gen- 
erated dynamically provided the three Sakharov condi- 
tions [i[ are fulfilled in the early universe: violation of 
baryon (or baryon minus lepton) number; violation of 
C and CP; and deviation from thermal equilibrium. In 
the standard model supplemented by heavy right-handed 
Majorana neutrinos, these conditions are naturally satis- 
fied for leptons. The Majorana mass term violates lepton 
number by two units. Complex Yukawa couplings of the 
right-handed neutrinos to leptons and the Higgs doublet 
induce CP violation. The rapid expansion of the universe 
causes a deviation from thermal equilibrium. Finally, the 
generated lepton asymmetry is converted to the observed 
baryon asymmetry by sphalerons @, @| . In other words 
the generation of the baryon asymmetry - baryogenesis 
- proceeds via the generation of a lepton asymmetry - 
leptogenesis 0. 

Many aspects of leptogenesis have been extensively in- 
vestigated. In particular, it has been shown that the 



* mathias.garny@ph.tum.de 



t andrcas.hohcncgger@mpi-hd.mpg.de 
t alcxander.kartavtsev@mpi-hd.mpg.de 
§ manfrcd . lindncr@mpi-hd . mpg.de 



CP-violating parameter and the efficiency of leptogene- 
sis are affected by the flavor structure of the neutrino 
Yukawa couplings [§4l6l|. In [13, [f8| it was demon- 
strated that the CP-violating parameter is resonantly 
enhanced if two of the heavy neutrinos have mass dif- 
ferences comparable to their decay widths. Medium ef- 
fects have also been addressed. In the hot and dense 
plasma the deviation of the CP-violating parameter and 
the thermal masses from their vacuum values plays an 
important role pjl, |2(j. In state-of-the-art calculations 
Boltzmann equations arc used to compute the asym- 
metry. Their applicability in the hot and expanding 
early universe can be checked using a first-principle ap- 
proach like the Schwinger-Kcldysh/Kadanoff-Baym for- 
malism. Some aspects of leptogenesis have been inves- 
tigated within this framework at different levels of ap- 
proximation in Minkowski space (2f| - [25| . These studies 
were motivated by the expectation that in the expand- 
ing universe filled with hot and dense plasma quantum 
effects, which are neglected in the canonical treatment, 
might play a crucial role. 

In this paper we investigate leptogenesis, and in par- 
ticular the vertex contribution to the CP-violating pa- 
rameter within the framework of nonequilibrium quan- 
tum field theory. Using the Kadanoff-Baym formal- 
ism as the starting point of our analysis, we derive 
quantum- corrected Boltzmann equations. We explicitly 
take medium corrections to the CP-violating parameter 
as well as the expanding background into account. The 
Kadanoff-Baym approach to the analysis of nonequilib- 
rium systems is technically considerably more involved 
than the canonical Boltzmann ansatz. For this reason, 
before applying it to realistic models of leptogenesis, here 
we study a simple toy model containing one complex and 
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two real scalar fields. It is defined by the Lagrangian 
C = hp^d^i - \Ml^ + d»bd„b - m 2 bb 

~ ^ ~ %^ hh ~ + C rest , * = 1, 2 , (1) 

where b denotes the complex conjugate of b. Despite 
its simplicity, the model incorporates all features rele- 
vant for leptogencsis. The real scalar fields imitate the 
(two lightest) heavy right-handed neutrinos, whereas the 
complex scalar field models the baryons. The U(l) sym- 
metry, which we use to define "baryon" number, is ex- 
plicitly broken by the presence of the last two terms, just 
as the B — L symmetry is explicitly broken by Majorana 
mass terms in phcnomcnological models. Thus the first 
Sakharov condition is fulfilled. The couplings gi model 
the complex Yukawa couplings of the right-handed neu- 
trinos to leptons and the Higgs. By rcphasing the com- 
plex scalar field at least one of the couplings <?i can be 
made real. If arg(<?i) ^ arg(<72) the other one remains 
complex and there is CP violation, as is required by the 
second Sakharov condition. In vacuum the vertex contri- 
bution to the CP-violating parameter is given by 




see AppendixlA"! Just as in realistic models, the required 
deviation from thermal equilibrium is caused by the rapid 
expansion of the universe. Thus the third Sakharov 
condition is fulfilled as well. Finally, the quartic self- 
interaction term in {TJ plays the role of the Yukawa and 
gauge interactions in established models - it brings the 
"baryons" to equilibrium. The renormalizability of the 
theory requires the presence of some additional terms, 
which are accounted for by £ re st- By appropriately 
choosing the corresponding coupling constants we can 
always make the contributions of these terms negligibly 
small. Since the physically interesting range for the gen- 
eration of the asymmetry is 0.1 Mi < T < lOMj, where 
Mi is the mass of the lightest heavy particle, the running 
effects cannot make these couplings large during the rel- 
evant period. 

Apart from the vertex contribution 0] to the CP- 
violating parameter discussed above, there is also a self- 
energy contribution [2til-[28|. In the Kadanoff-Baym for- 
malism the analysis of the former is rather independent 
from the analysis of the latter. For this reason, in this 
paper, we consider only the vertex contribution, whereas 
the self-energy contribution will be addressed in [29[ . To 
make the discussion more transparent we give the techni- 
cal details in the appendices, whereas in the main body of 
the paper we discuss qualitative features of the employed 
approach and present the results. 

(i) As we argue in Sec.|Tll the formalism is free of the 
double-counting problem typical for the canonical 
Boltzmann approach. In other words the struc- 
ture of the equations automatically ensures that 



the asymmetry vanishes in thermal equilibrium and 
no need for the real intermediate state subtraction 
arises. 

(ii) Our result for the vertex contribution to the CP- 
violating parameter, presented in Sec. Mil differs 
from that obtained in the framework of equilib- 
rium thermal field theory by replacing the zero- 
temperature propagators with finite temperature 
propagators in the matrix elements of the Boltz- 
mann equation fl9l . I20I ] - the medium corrections 
are only linear in the particle number densities. 
For scalars the medium effects always increase the 
CP-violating parameter, which in turn leads to an 
enhancement of the generated asymmetry. 

(iii) By comparing the CP-violating parameters ob- 
tained by using the Maxwcll-Boltzmann (MB) and 
Bose-Einstein (BE) statistics, we find that quan- 
tum statistical effects play a considerable role. As 
we argue in Sec. lIIIi the medium effects increase 
the CP-violating parameter by a factor of at most 
two in the Maxwell-Boltzmann approximation. At 
high temperatures, the increase is up to an order of 
magnitude larger when Bose enhancement is taken 
into account. 

In Sec. lIVI we present numerical solutions of the quantum- 
corrected Boltzmann equations, and discuss the quanti- 
tative impact of medium effects on the final asymmetry 
within the toy model. Finally, in Sec.|Vl we summarize 
our results and present our conclusions. 

II. NONEQUILIBRIUM DYNAMICS 

To calculate the asymmetry generated at the epoch of 
leptogenesis one usually employs generalized Boltzmann 
equations for the one-particle distribution functions of 
the different particle species [l|]: 

P a v a f 4 ,(x,p) = i / dn*»dng . . . dnfdnf . . . 

x (2ir) 4 5(p + p a +Pb.-.—Pi~Pj) 
x [\M\ 2 i+j+ ^ +a+b Jdi ■ ■ • (1 ± fa)(l ± Ml ± U) 
f a f b U...(l±f i )(l±f j )..}. (3) 

The fourth line in Eq. (|3| describes the decrease in num- 
ber of species ip due to the scattering (or decay) process 
ip + a + b . . . — > i + j + ... and is usually referred to as 
the loss term. The third line describes the increase in the 
number of tp due the process + — > ijj + a + b . . . and 
is referred to as the gain term. The Dirac 5 function in 
the second line enforces energy-momentum conservation 
in each individual process, whereas the invariant phase- 
space elements dH 3 ensure that the resulting expression is 
a Lorentz scalar. The probabilities of the decay and scat- 
tering processes are usually calculated in vacuum, which 
is inconsistent with the nonzero particle number densi- 
ties. Moreover the canonical approach is plagued with 
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the double- counting problem. For instance, in the canon- 
ical approach, the scattering process £h ipi £h is 
equivalent to the inverse decay (£h, £h—> of the heavy 
Majorana neutrino followed by the decay (tpi —> £h,£h) 
if the intermediate heavy neutrino tpi is on-shell. That 
is, the same contribution is counted twice. As a conse- 
quence, a nonzero asymmetry is generated even in ther- 
mal equilibrium. The problem is accounted for by the 
real intermediate state subtraction procedure. Since the 
scattering amplitude is calculated in vacuum, one cannot 
assign a distribution function to the heavy neutrino ipi. 
For this reason the resulting collision terms are difficult to 
interpret. Only by assuming Maxwell-Boltzmann statis- 
tics and by integrating the resulting Boltzmann equations 
can one derive rate equations that manifestly lead to zero 
asymmetry in equilibrium. 1 

Boltzmann equations rely on the concept of on-shell 
particles with constant mass. The spectral function 
G p (X,p) of a particle whose mass does not change as 
it propagates along a geodesic is orthogonal to its four- 
momentum [30l ]: 

p a V a G p (X lP ) = 0. (4) 

For a pointlike on-shell particle the spectral function is 
zero off-shell and diverges on-shell: 

G p (X,p) = 2 7 rsign( Po ) ( 5(p> Q - M 2 ) . (5) 

Instead of using one-particle distribution functions, as 
the quantities describing the statistical properties of the 
system, we can use the statistical propagator Gf(X,p): 

G F (X,p) = [U(X,p) + \] G p (X,p) . (6) 

Written in terms of G p (X,p) and Gf(X,p), a Boltzmann 
equation takes the form 

p a V a G F (X,p) 

= l[U < (X lP )G > (X,p) - G<(X,p)n > (X,p)] , (7) 

where = Gf ± \ G p . Comparing Eqs. and © we 
see that the quantities correspond to the loss and 
gain terms. 

The transformations made so far simply seem to 
amount to a change of notation. However, the situa- 
tion changes when one realizes that Eqs. (j4]) and ([7|) 
coincide with the equations that can be derived in a cer- 
tain approximation from the system of Kadanoff-Baym 
equations and that the quantities can be identified 
with the self-energies. Despite the close similarity there 



Strictly speaking, the system Q is incorrect in general if it in- 
volves unstable particles, because the matrix elements, computed 
in perturbation theory, do not necessarily meet the requirements 
for transition amplitudes in systems of Boltzmann equations. In 
cases such as the present (leptogenesis) subtle modifications (RIS 
subtraction) are necessary to obtain consistent results. 



is an important difference: the self-energies calculated 
in the framework of the Kadanoff-Baym formalism differ 
from the gain and loss terms obtained in the canonical ap- 
proach. As we will show in the following, the self-energies 
consistently take medium effects into account and the re- 
sulting equations are free of the double-counting problem. 

The Kadanoff-Baym formalism may be viewed as top- 
down approach: starting from the complete evolution 
equations for the two-point functions, it is possible to 
derive kinetic equations in a systematic wa y by a pply- 
ing a number of well-known approximations [3 Xl — 1341 ] . We 
will refer to these equations as quantum- corrected Boltz- 
mann equations. Within the Kadanoff-Baym formalism, 
both the overall structure of the equations as well as the 
scattering and decay rates are derived sclf-consistently 
from a common starting point based on the in-in or 
Schwinger-Keldysh description of nonequilibrium quan- 
tum fields. In contrast to that, within the canonical 
bottom-up approach, the scattering and decay rates are 
extracted from the S-matrix (i.e. based on the in-out for- 
malism), and are then inserted into Boltzmann equations. 
We expect that both approaches are equivalent in the 
zero-temperature/zero-density and small coupling limit, 
where the mean-free path is large compared to the mi- 
croscopic interaction length scales. In this regime, it may 
be considered that the top-down approach discussed here 
adds support to the usual bottom-up formalism. Addi- 
tionally, as mentioned above, the top-down approach is 
free of the double-counting problem and allows one to 
explore the medium effects within nonequilibrium quan- 
tum field theory. This is especially relevant for the CP- 
violating terms, since these contain loop diagrams. 



A. Kadanoff-Baym equations 

The system of Kadanoff-Baym equations for the spec- 
tral function and the statistical propagator is usually for- 
mulated in coordinate space. For the complex scalar field 
b, which corresponds to baryons in our model they read 
(see Appendix lBj) 

y° 

[n x + m 2 (x)]D F (x,y) = J S i zS F (x,z)D p {z,y) 



x° 

- J@ 4 zZ p (x,z)D F (z,y), (8a) 


v° 

[U x +m 2 (x)]D p {x,y) = J ^z^ p (x, z)D p (z,y) . (8b) 

x° 

As is clear from the terminology, the spectral function 
D p contains information about spectral properties of the 
system, whereas the statistical propagator Dp contains 
information about its state. The spectral and statistical 
self-energies, E p and Ej?, as well as the effective mass 
m 2 {x) = m 2 + Y,i oc (x), which contains the local self- 
energy, carry information about the interactions in the 
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system. They describe scattering and mean-field effects, 
respectively. The invariant volume clement, 

= ^-gd 4 z , g = dct g^ v , 

ensures that the Kadanoff-Baym equations (JSJ) can be 
applied to the analysis of out-of-cquilibrium dynamics 
not only in Minkowski, but also in a general curved space- 
time. Particularly interesting for our purpose is the case 
of the expanding early universe J35| ■ 

Let us list some of the qualitatively important fea- 
tures. First of all, Eqs. (j8| are written in terms of re- 
summed propagators, i.e. they take into account the full 
series of daisy and ladder diagrams (see e.g. [HI). Sec- 
ond, the characteristic memory integrals on the right- 
hand side integrate over the full time history of the evo- 
lution. In other words the Kadanoff-Baym equations are 
non-Markovian, i.e. are not local in time. It is very im- 
portant that the Kadanoff-Baym equations do not rely 
on the concept of quasiparticles and their collisions in 
the plasma. In other words, they are free of any possi- 
ble uncertainties associated with definition of quasiparti- 
cle excitations in the hot plasma of the rapidly expand- 
ing universe. This property makes the Kadanoff-Baym 
equations a prime candidate for the analysis of lepto- 
genesis. If the quasiparticle picture is applicable, they 
account for the time-dependence of the quasiparticle pa- 
rameters. In particular, an effective time-dependent mass 
and width induced by the interactions of the system can 
be extracted from the Wigner-transform of the spectral 
function [33] ■ Finally, for weakly coupled systems close to 
thermal equilibrium the Kadanoff-Baym equations can 
be reduced to the Boltzmann equation [3ll434| , which we 
have briefly discussed above. 

The Kadanoff-Baym equations for a system of n real 
scalar fields read [H, US] 

y° 

[□, + Mf]G^(x, y) = f® 4 z V*{x, z)G k /(z, y) 
o 

- X f®*zILf(x,z)G k J(z,y), (9a) 
o 

y° 

[□, + MfC^U. y) = f&zl$(x, z)G k /(z, y), (9b) 

x° 

where the propagators and self-energies are now n-by-n 
matrices. The off-diagonal components of the propaga- 
tors and self-energies describe the mixing of the fields. 

Equations (|8]) and ([9]) are valid for a system of one 
complex and n real scalar fields with arbitrary inter- 
actions. Here, we consider n = 2. The information 
about the particular form of the interactions is encoded 
in the corresponding self-energies S(x,?/) and Ii l3 {x,y). 
The latter ones can be derived by functional differentia- 
tion of the two-particle-irreducible (2PI) effective action, 
T2 [G, D] , with respect to the two-point correlation func- 
tions D(y,x) and G 3l {y,x), see [301 an d AppendixlBl for 
more details. The effective action is given by the sum of 



all 2PI diagrams with vertices as given by the interaction 
Lagrangian and internal lines representing the complete 
connected two-point functions [38|. Of course the num- 
ber of 2PI diagrams contributing to the effective action 
is infinite in any theory and the infinite sum must be 
truncated in a suitable way. To obtain a qualitative sim- 
ilarity between the toy model and the established models, 
we must take into account processes which generate and 
washout the asymmetry. The minimal set of relevant 2PI 
contributions is presented in Fig.[TJ Clearly, after cutting 




(a) (b) (c) (d) 



FIG. 1: Relevant two- and three-loop contributions to the 
2PI effective action. 

any two lines of the diagrams (a) — (d) they still remain 
connected. To understand which physical processes are 
described by the above diagrams in the Boltzmann ap- 
proximation one has to discriminate between local and 
nonlocal contributions. The local diagram (a) generates 




FIG. 2: Processes described by the 2PI diagrams in Fig.Q] 

the mean-field correction E loc (a;) to the effective mass 
m 2 (x) of the field: 

E loc (x) = XD(x,x) . (10) 

The nonlocal diagrams (b)-(d) describe scattering and 
decay processes, which can be identified by cutting the 
diagrams into two pieces by drawing a connected line in 
all possible ways as indicated in Fig. [5] (see also [39|). 
Cutting diagram (b) , we obtain squares of tree- level am- 
plitudes of bb — > bb and bb — > bb scattering processes. 
Analogously, cutting diagram (c), we get squares of the 
tree-level amplitudes of (/', — > bb and (/', — > bb decay pro- 
cesses. In the canonical approach the decays of the heavy 
real scalars are CP-conserving at tree-level. To leading 
order the vertex CP-violating parameter e is generated 
by interference of the tree-level and one-loop vertex de- 
cay amplitudes. In the Kadanoff-Baym formalism the 
leading vertex CP-violating contribution is described by 
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FIG. 3: Interference of the tree-level and one-loop vertex de- 
cay amplitudes. 



diagram (d) . Cutting it in the way presented in Fig. [31 we 
obtain the product of the tree-level and one-loop vertex 
amplitudes. 2 

Instead of calculating the spectral and statistical com- 
ponents of the self-energies, it is easier to calculate the 
Wightman components £^(x, y) = Ep(:r, y) T |S p (x, y). 
These can be identified with the gain and loss terms, as 
mentioned earlier in this section. After some simple but 
tedious algebra (see AppcndixICl for more details), we 
obtain the self-energies corresponding to the diagrams in 
Fig.ffl 

Ef(x,y) = -^X 2 D k (x,y)D^(x,y)D^(y,x), (11a) 
SjW) = -<H9jG%(v,*)D$(y t x), (lib) 
Xf(x,y) = -g i9j g* m g* n f @*v^u 
[ D R (y, v)D F (u, »)%(«, x)G%{y, u)Gf l (x, v) 
+D R (y, v)D A (u, v)D^(u, x)G%(y, u)G^ n (x, v) 
+D F (y, v)D R (u, v)D$(u, x)G%(y, u)G^ n (x, v) 
+%(y, v)D F (u, v)D A (u, x)G%(y, u)GT{x, v) 
+%(y, v)D A (u, v)D F (u, x)Gl(y, u)GT(x, v) 
+D $ (y, v)D R (u, v)D A (u, x)G^(y, u)G^ n {x, v) 
+D R (y, v)D^ («, v)D A (u, x)G| (y, u)G^ n (x, v) 
+D^(y, v)D$(u, v)D^{u, x)G%{y, u)G% n (x, v)] . (11c) 

Diagrams (6) and (c) in Fig. [1] induce contributions T,^ 
and £( c ) which contain only the Wightman two-point 

correlation functions D> and G l -i . As we will show later, 

< < 

in the Boltzmann approximation these correspond to the 
on-shcll initial and final states. To write the last term in 



compact form, we have introduced the retarded and ad- 
vanced propagators, D R (x,y) = 9(x° — y°)D p (x,y) and 
D A (x,y) = ~0{y° — x°)D p (x,y), respectively. The first 
six terms 3 of describe the one-loop correction to 
the decay width. The combinations of the statistical, re- 
tarded and advanced propagators in Eq. (jllcp correspond 
to the three internal lines in the loop, whereas the ^ com- 
ponents again correspond to the on-shcll initial and final 
states. 

Since we do not consider quartic interactions of the 
real scalar fields, there are no local corrections to their 
masses. It is for this reason that Eqs. Q contain only 
the bare masses Mf of the fields. The first nonlocal term 
IT( C ) describes the decay of the heavy scalar at tree-level: 



4 c)u (x,y) = -\ gi g*D\(x,y) - y* 9j D^(y,x) , (12a) 

n[f J (z,y) = -fagjg^* J 0Wu 

[ GT L (v,u)D^ (a, v)D k (x, u)D R (y, v)D R (y, u) 
+GT (v,u)D^ (x, v)D^ (x, u)D R (y, v)D F (y, u) 
+GT(v,u)D% (x, v)D^ (x, u)D F (y, v)D R (y, u) 
+G% n (v,u)D R (x, v)D R (x, u)D $ (y 7 v)D $ (y, u) 
+GT(v,u)D R (x, v)D F (x, u)D $ (y, v)D^(y, u) 
+G A n {v,u)D F (x 7 v)D R (x, u)D $ (y, v)D $ (y, u) 
+G f l (v,u)D^ (x, v)D R (x, u)D R (y, v)D $ (y, u) 
+G^ n (v,u)D R (x, v)D^ (x, u)D$(y, v)D R (y, «)] 

[ G% n (v,u)D $ (v, x)D $ (u, x)D A (v, y)D A (u, y) 
+ ... ]. (12b) 

The first six terms 4 in each of the two square brackets 
of ir d ' describe the one-loop corrections to the scatter- 
ing width. Their structure is very similar to that of the 
first six terms of Eq. (|llcj) and the combinations of the 
statistical, retarded and advanced propagators again cor- 
respond to the three internal lines in the loop. 

The Kadanoff-Baym equations ([5]) and © together 
with the expression for the self-energies (fTTj) and (fT2"]) 
form a closed system of integro-diffcrcntial equations. Its 
solutions carry full information about the spectral and 
statistical properties of the system, including information 
about the generated asymmetry at each instant of time. 



2 In addition, there are two other ways to cut this diagram, which 
are denoted by the ellipses in Fig. [2] They describe scattering 
processes bb — ¥ bb, bb — > ipi^Pj an d ipib — > ipjb. Note that the 
three-loop 2PI diagram only describes the interference of the s, t 
and M-channel scattering amplitudes: M a t, M su , Mtu- The miss- 
ing topologies, which generate the M aa , Mtt, and M uu terms, 
appear only upon use of the extended quasiparticle approxima- 
tion. This analysis is beyond the scope of this paper and will be 
presented in |29l . 



3 The seventh term in (II let describes the scattering process tpib —¥ 
ipib. This is clear from the fact that it contains one and 
two Gjj two-point functions, i.e. one "external" complex scalar 
and two "external" real scalars. Similarly, the eighth term of 
Eq. Ullcjl describes the scattering process bb — > bb, because it con- 
tains three two-point functions, i.e. three "external" complex 
scalars. 

4 The seventh and eighth terms of Eq. 112bt describe the scattering 
processes i/^fc — > ipjb and bb —¥ ipi^Pj- 
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B. Quantum-corrected Boltzmann equations 

Despite all advantages, the full Kadanoff-Baym equa- 
tions are relatively rarely used for the analysis of out- 
of-equilibrium processes, partially because of the com- 
plexity of the numerical solution. 5 In this section we de- 
rive kinetic equations starting from the above Kadanoff- 
Baym equations by applying agradient expansion and a 
quasiparticle approximation |3ll434l . l39l . |49| . [Hoj . The re- 
sulting quantum-corrected Boltzmann equations can be 
directly compared to the canonical equations and are eas- 
ier to solve numerically. 

The quantum-corrected Boltzmann equations are ap- 
plicable to weakly coupled systems of (quasi) particles 
that have a width which is small compared to their mass 
and that evolves slowly compared to the microscopic in- 
teraction time scales. We expect that these conditions 
are satisfied for thermal leptogenesis, where the devi- 
ations from equilibrium arc moderate in general. In 
thermal equilibrium the two-point correlation functions 
D(x,y) and G t: '(x,y) depend only on the relative coor- 
dinate, s = x — y, and are independent of the center 
coordinate 6 , X = ^(x + y). Having these equilibrium 
considerations in mind, we trade the variables x and y 
for the new arguments X and s: D(x,y) — > D(X,s). 
Out of equilibrium the two-point functions depend on 
both the relative and center coordinate. If, however, the 
deviation from equilibrium is small, one can perform a 
gradient expansion of the correlation functions and the 
self-energies in the vicinity of X keeping only the leading 
terms. This results in a system of equations describ- 
ing the slow, AT-dcpcndent dynamics of the statistical 
properties of the system, see AppcndixlBl for more de- 
tails. Performing the gradient expansion, which is the 
first step in the derivation of the Boltzmann equation, 
we replace the non-Markovian evolution equations by a 
system of Markovian ones. Therefore, when truncating 
at first order in the gradient expansion, we neglect the 
memory effects. The fast dynamics associated with the 
relative coordinate is responsible for the spectral prop- 
erties of the system, which are conveniently described in 
the momentum representation. Performing the Wigncr 
transformation (see Appcndix lB]) . we trade the relative 
coordinate s for a coordinate p in momentum space: 
D(X,s) — > D(X,p). To perform the Wigner transfor- 
mation we have to send the initial time to minus infinity 
which means that we neglect the effects of initial corre- 
lations. The next step in the derivation of Boltzmann 
equations is the Kadanoff-Baym ansatz. It relates the 



statistical propagator and the spectral function by a gen- 
eralized fluctuation-dissipation relation, 



5 See e.g. HHH]. 

6 The above definitions of the relative and center coordinates are 
valid only in Minkowski space-time. In a general space-time the 
center coordinate is defined as coordinates of the center (fx ) of 
the geodesic connecting x and y, X a = £ a (fx), whereas the 
relative coordinate is proportional to the length of the geodesic 
between the two points, s a = ti a (?x)(?;c ~ f») l5ll . 



D F (X,p)= [f{X,p) + \]D p {X,p), 



(13) 



where f{X,p) is the one-particle distribution function. In 
the equilibrium limit, the Kubo-Martin-Schwinger con- 
dition [s^, ensures that f(X,p) converges towards a 
Bose-Einstein distribution function. The final step is the 
quasiparticle approximation, where one replaces the ex- 
act smooth spectral function by a Dirac 6 function peaked 
on the mass-shell of the quasiparticles. The resulting 
Boltzmann-like equation, which describes the time evo- 
lution of the particle distribution function, reads 

[p a V a f(X lP )}D p (X,p) 

= |[D < (X,jj)E > (Jf,p) - D > (X,p)X < (X,p)} . (14) 

The analogous equation for antiparticlcs is given by 

[p a V a f(X lP )}D p (X,p) 

= p)£>(X,p) - D>(X,p)£<(X,p)] , (15) 

where D^(X,p) = t^(X,p) = %(X,-p). 

Note again that by making the approximations which 
have led to (fT4"|) and (fTBj) we have neglected the memory 
effects and the effects of the initial correlations. As a 
result, the quantum-corrected Boltzmann equations are 
Markovian, i.e. local in time. 

To obtain a closed system of quantum-corrected Boltz- 
mann equations, we must also Wigner-transform the self- 
energies (fTTj) (sec AppcndixiDl for more details). The lat- 
ter encode the scattering and decay rates including quan- 
tum noncquilibrium effects. By employing the relations 

D < {X,p) = f{X,p)D p (X,p), 
D > (X,p) = [l + f(X,p)}D p (X,p), 

which follow directly from Eq. (|13|). we can then rewrite 
Eqs. (|14I15[) in a way resembling the usual form of Boltz- 
mann equations, including the correct quantum statisti- 
cal factors. 

Within the 2PI three-loop approximation, we find that 
there are two physically distinct contributions to the self- 
energy. The first one, corresponding to H^ b \ describes 
CP-conserving two body scatterings, bb — > bb, at tree- 
level: 

s|(x,p) = -±a 2 / <m pl <m P2 <m P3 {2TT)H3{p+p l 

-P2 - p 3 )% (X, Pl )D^ (X, p 2 )D^ (X, p 3 ) , (16) 

where the invariant volume element in momentum space 
is defined by 



dll 



1 



d 4 p 



(17) 



and 5 9 (p) e 
the S function [3C 



~9x (2^) 4 ' 
8{p) is the covariant generalization of 
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The second contribution, given by the sum of £( c ) and 
describes decay processes ipi — s- bb and ipi — > bb at 
tree- and one-loop level: 

£|(X,p) = -|. gi | 2 /dn pi dn P2 (27r) 4 ^(pi - pa -p) 

x [l + A|(X,p 1 ,p 2 )]G|(X,p 1 )%(X,p 2 ). (18) 

The newly introduced function A^(X,pi,p2) takes into 
account the one-loop corrections to the decay width and 
is given by 

Al(X,p uP2 ) = \ 9j \ 2 ( 9 4A f dU kl dU k2 dU k3 
\ 9i 9j / J 

x (27r) 4 <F(pi + h + k 2 ){2ll)H 9 {k 2 - k 3 + P2 ) 

+ D A (X,k 1 )D R (X,k 2 )G^(X,k 3 ) 

+ D F (X,k 1 )D A (X,k 2 )G J i(X,k 3 )]+ c.c. (19) 

Proceeding in the same way, we derive quantum- 
corrected Boltzmann equations for the distribution func- 
tions of the real scalar fields, which is a two-by-two dif- 
ferential matrix equation. The off-diagonal components 
of the correlation functions are generated dynamically by 
the exchange of two complex scalars and are therefore of 
order g 2 . The one-loop vertex terms, which generate the 
CP-violating parameter, are proportional to the fourth 
power of the coupling constant. Therefore the contribu- 
tion of the off-diagonal terms to the vertex CP-violating 
parameter is of the order of g 6 . Here, we limit ourselves 
to terms of at most fourth power in the coupling constant 
and therefore we can neglect the off-diagonal terms in the 
corresponding matrix equation. The resulting equations 
coincide then with those derived in (30j : 

[p a v a f^(x : p)}G i ;(x : p) 

= i[G»(X,p)U«(X,p) - G"(X,p)U«{X, P )] . (20) 

Note that we have in fact used this diagonal approxima- 
tion in Eqs. (fT%|) and (flT?)) . The Wigncr-transform of the 
self-energy (|12j) is given in the same approximation by 

uf\x, P ) = -i|.g 4 | 2 /dn J)1 dn P2 (2^) 4 ^( Pl + P2 - P ) 

x {[1 + Atp(X,p,p 2 )] D^(X, Pl )D^(X, P2 ) 

+ [l + ^{X,p,p 2 )}D^{X,p x )D^{X,p 2 )}. (21) 

The second line of (|2T|) describes the process tp — > bb. 
The one-loop correction to this process is given by 

&%{x,p,p 2 ) = \ 9l \ 2 ( g 4A I dn kl du k2 du k3 

\ 9i 9j / J 

x (2tt) 4 5S{ p + kt + k 2 ){2n) 4 S 9 {k 2 - k 3 + p 2 ) 
x [D R (X,k 1 )D R (X,k 2 )G^{X,k 3 ) 

+ D R {XM)DF{X,k 2 )G 3 i{X,k 3 ) 

+ D F {X,k 1 )D R {X,k 2 )G%(X,k 3 )]+c.c. (22) 



The third line of (|2"Tj) describes ip —> bb process and the 
corresponding one-loop contribution is related to (|22[) by 
A;(X,p 1 ,p 2 )^A;(X,- Pl ,-p 2 ). 

A very important feature of the expressions for the self- 
energies, Eqs. (TTSJ) and (|2"Tj) . is that the loop corrections 
A£ and A^, appear as overall factors on the right-hand 
sides of the corresponding quantum-corrected Boltzmann 
equations. Therefore, when using the conventional ap- 
proximations and integrating Eqs. (|14ll5p . we obtain the 
structure 7 



d t (n b , n- b ) oc (1 ± e,)(n. t - nf) .... (23) 

To obtain an equivalent result in the canonical approach 
one explicitly needs to apply the MS subtraction proce- 
dure. This means that, here, the structure of the equa- 
tions automatically ensures that no asymmetry is gen- 
erated in thermal equilibrium. Stated differently, the 
Kadanoff-Baym formalism is free of the double- counting 
problem and no need for MS subtraction arises. 

In the homogeneous and isotropic early universe the 
canonical Boltzmann equations conserve the linear com- 
bination 2n^ 4 + rib + n>i °f particle numbers in a comov- 
ing volume. However, the conservation of this quan- 
tity is accidental, i.e. not guaranteed by a symmetry 
of the underlying Lagrangian. It is not conserved by 
the full Kadanoff-Baym equations 8 (see [47[ for another 
example). This is also true for the quantum- corrected 
Boltzmann equations which we have derived from the 
Kadanoff-Baym equations. To see this one should add 
Eqs. (TTJ} and (JTSJ) as well as Eq. (|2DJ) multiplied by two 
and use the explicit expressions for the self-energies (fT8)) 
and (|2Tj) . Although the expressions for the loop correc- 
tions ([TO)) and (f2"2"|) are similar, they are not equal. This 
results in a small time-dependence of the quantity men- 
tioned above (see also Appendix lG|). 

The system of Boltzmann equations (TH|) , (JT5)) and (f2"Tj)) 
together with the Wigner-transforms of the self-energies 
(fT6|l . (fl8|) and (f2Tj) form a closed system of differential 
equations that can be solved numerically. The solutions 
describe the phase-space distributions of quasiparticle ex- 
citations at each instant of time t = X°. 



7 For the numerical analysis we use the full Boltzmann equa- 
tion, since the approximations required to obtain integrated 
rate equations are not appropriate within the toy model (see 
Scc. lIVII . Note, however, that the latter in fact contain averaged 
effective CP- violating parameters (see Sec. lUII and AppcndixlFl . 
£i -> (ei)- 

8 ft is important to note that the 2PI approximation scheme guar- 
antees that the total energy-momentum tensor (including "po- 
tential" energy due to interactions in the system) is covariantly 
conserved |33l. l54j. In other words, the nonconservation of the 
total number of particles does not violate the fundamental con- 
servation laws. 
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III. CP- VIOLATING PARAMETER 

Comparing the Boltzmann equations for particles and 
antiparticles, (fFl)) and (fT5"]) , we see that the dynami- 
cal generation of the CP asymmetry is only possible if 
^(X,p) ± %{X,p). Since %{X,p) = %(X,-p), in 
the diagonal approximation, this is equivalent to the re- 
quirement that A\{X,p\,p 2 ) 7^ Al(X,—pi,—p 2 ). The 
CP-violating parameter can then be defined as 

ei (X, pi,p 2 ) = l [Aj (X, Pl ,p 2 )-Ai (X, - Pl , -pa)]. (24) 

Using properties of the Wigner transforms of the statis- 
tical, retarded and advanced propagators, we find that in 
a medium that is (approximately) baryosymmetric 9 it is 
given by 



ei(X, Pl ,p 2 ) = |.9j| 2 Im 



9,9, 



dn kl dn k2 diL k3 



,9193. 

x {2^)H a { Pl + fci + fc 2 )(2^) 4 ^(fc 2 - k 3 +p 2 ) 

x [D P (X, k x )D F {X, k 2 )G{ j (X, k 3 ) + {k! o k 2 } 

+ D h (X, ki)D P (X, k 2 )Gf(X, k 3 ) + {h o k 2 } 

+ D p (X,k 1 )D h {X,k 2 )G j J{X,k s ) - {h o fc 2 }] , (25) 

where D h (X,p) = RcD R (X,p) = RcD A (X,p) and 
is defined analogously. The quasiparticle approximation 
together with the Kadanoff-Baym ansatz enforces the 
spectral functions and the statistical propagators to be 
on mass-shell. On the contrary, the real parts of the 
retarded propagators, Dh and Gh, vanish on-shell (see 
AppcndixlBl for more detail). In other words, in the pro- 






FIG. 4: Graphical representation of the terms in Eq. (I25|) . 

cesses contributing to the CP-violating parameter, two 
of the intermediate lines are on-shell and one line is off- 
shell; this is shown in Fig. [H It is also interesting to 
note that in each term on the right-hand side of (|25|) 
only one of the internal lines is "thermal", i.e. explic- 
itly depends on the one-particle distribution function. In 



other words the medium corrections are only linear in 
the particle number densities. One could come to this 
conclusion even without explicitly calculating the CP- 
violating parameter. The self-energy (|llc| is a product 
of five two-point correlation functions. As far as the de- 
cay is concerned, two of them describe on-shell "exter- 
nal" states. The two integrals over the space-time coor- 
dinates, u and v, turn two other functions into retarded 
and (or) advanced propagators, which do not explicitly 
depend on the particle number densities. The remaining 
function turns out to be given by the corresponding sta- 
tistical propagator and explicitly depends on the particle 
number density. 

Typically, one is interested only in the total generated 
asymmetry and solves rate equations for the total particle 
number densities. They are obtained by integrating the 
left- and right-hand sides of the Boltzmann equations 
([14"]) and (TToj) over phase space. In the corresponding 
gain and loss terms, 

/ dn p dll pi dTL P2 (2tt) 4 ^ (pj - p 2 - p) [1 ± e t (X, pi , p 2 )] 
x G|(A,p!)%(A,p 2 )%(A,p) , (26) 

we can perform the transformation p O p 2 and take 
a sum of the initial and final expressions, so that in 
Eq.®: 

ei(X,p x , P2 ) -> \{e i {X,p 1 ,p 2 ) + e i {X,p 1 ,p)]. (27) 

In ei(X, P i,p), we transform the variables k\ -R- k 2 
and k 3 — > — fc 3 in addition. The spectral function 
of a real scalar field is antisymmetric: G % HX, — k\) = 
—Gp(X,ki). Collecting all the terms, we find that due 
to the antisymmetry only the first two terms (the third 
line) in (f2"5j) contribute to the right-hand side of (|27|). 
whereas the other terms cancel out. In other words, only 
the first two terms of (f2"5j) contribute to the total CP- 
asymmetry, whereas the other four terms do not. An 
explicit calculation shows that (at least in the homoge- 
neous and isotropic universe) these four terms also do 
not contribute to the gain and loss terms. Diagrammati- 
cally this means that only decays followed by a scattering 
contribute to the CP-violating parameter, see Fig.H^,. 

As we have already mentioned, it is important that 
the CP-violating parameters are identical 10 for the 
gain and the loss terms: this means that the structure 
of the quantum- corrected Boltzmann equation automat- 
ically ensures that the asymmetry vanishes in thermal 
equilibrium. Let us also note that there is a clear dis- 
tinction between the initial, final and intermediate states: 
the former ones are described by the Wightman functions 



We will see later that the CP-violating parameter, defined in 
this way, is different for particles and antiparticles if the corre- 
sponding distribution functions are different. Since the expected 
asymmetry is small, this is only a second order effect and can 
be neglected in the present work. The condition of almost zero 
asymmetry is certainly satisfied if the CP- violating parameter is 
calculated in vacuum, as it is the case in the canonical approach. 



It is important to note that the transformation properties under 
C , P and T of the self-cncrgics obtained from the Schwingcr- 
Keldysh/Kadanoff-Baym formalism cannot be identified with 
those of the S-Matrix elements appearing in the canonical in- 
out formalism. 
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and G^, whereas the latter ones are described by the 
retarded, advanced and (or) statistical components of the 
two-point functions. 

Applying the quasiparticlc approximation and the 
Kadanoff-Baym ansatz in Eq. (|2"5|) , we obtain the follow- 
ing expression for the CP-violating parameter which is 
one of our central results: 



ei(pi,P2) = - 



1 



1.9,1 



-Im 



9r9n 



Sir Mf \gr 9jj 

dn i + J(E kl ) + f(E k2 ) 



4tt Mf/Mf + + 



cos ( 



(28) 



where Pfc 12 are the energies of the intermediate toy bary- 
ons as a function of pi, p 2 and the angle variables, and 
we have omitted the time-space coordinate X to shorten 
the notation. The CP-violating parameter is a sum of 
vacuum and medium contributions. Integrating the vac- 
uum contribution over the solid angle, we obtain the 
standard expression for the CP-violating parameter, see 
Eq. ([5]). The thermal contributions are proportional to 
the one-particle distribution function, which is positive. 
Hence, for scalars the CP-violating parameter is always 
enhanced by the medium effects. 

Because of the fact that the intermediate toy baryons 
propagate with respect to the rest frame of the thermal 
bath, the CP-violating parameter depends in each indi- 
vidual decay on the phase-space distribution of the de- 
caying particle and the decay products. Using results of 
AppcndixlEl we obtain for the energy of the intermediate 
complex scalars: 



E, 



fel,: 



h [Ei + |pi | (sin 8 cos ip cos 5' + cos 9 sin 5')} , (29) 



where 9 and ip are elements of the solid angle O and 
the angle 8' depends on momenta of the initial and final 
states: sinS' = (|p 3 | - |p 2 |)/|pi|. 

In the limit of almost equal one-particle distribution 
functions of particles and antiparticles, / and /, the CP- 
violating parts of A^, do not contribute to the Boltzmann 
equations for the real scalars, just as in the canonical 
approach. For this reason, we do not consider it here. 

In order to estimate the size of the medium corrections, 
we consider the hierarchical limit of the heavy scalar mass 
spectrum, Mi -C M 2 . As in standard leptogenesis, we as- 
sume that the asymmetry is predominantly generated by 
the decay of the lighter scalar. By expanding Eq. (f2"5|) in 
Mf/Mf, we obtain a simplified expression for the rele- 
vant CP-violating parameter ex, 



ei(Pi)P2) 



[1 + JiF {f(E kl ) + f(E k2 )}] , (30) 



where e\ ac is the CP-violating parameter in vacuum 
given in Eq. (J2J). Exploiting the k\ O k- 2 symmetry and 
integrating over the full solid angle we find that the asym- 
metry depends on the absolute value of pi only. That is 
ei(Pi,P2) = ex(|Pi|), where 



ei(|p|)=er 



1 



E max /2 



r\P\ JE mm /2 



f(E)dE 



(31) 



and E max = E x + r|p| and E mm = E x - r|p| are the 
largest and smallest kinematically allowed energies of 
the light scalars produced in the decay ipi — s- bb. Here 
Ei = (Mi +p 2 ) 5 and |p| denote the energy and momen- 
tum of ipi in the rest frame of the medium, respectively. 
Furthermore, we have also included a nonzero "baryon" 
mass m for completeness, which enters via the parameter 
r= (1 -4m 2 /Mi 2 )s. 
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FIG. 5: Effective CP-violating parameter ei(|p|) in medium 
obtained from the Kadanoff-Baym formalism. The shaded 
areas correspond to the range 0.25 < |p|/T < 4 of momenta 
|p| of the decaying particle ipi — > bb/bb with respect to the 
rest frame of the medium. Here we assumed a thermal Bose- 
Einstein (BE) and Maxwell-Boltzmann (MB) distribution for 
6/6 with vanishing chemical potential for illustration. In the 
low-temperature limit (NR), the vacuum value is approached. 
In the high-temperature limit (UR), the CP- violating param- 
eter is enhanced for bosons. We also show the thermally aver- 
aged CP-violating parameter (ei) for the BE (red long-dashed 
line) and MB (blue dashed line) cases, as well as the result 
that would be obtained in thermal field theory (green dotted 
line). 

The medium correction depends on the one-particle 
distribution function of the light scalars and the masses 
of the particles. As expected they vanish if / = 0. We 
emphasize that the upper expression is valid even if the 
light scalars were out of equilibrium. Nevertheless, since 
we expect the light scalars to be close to kinetic equilib- 
rium at all times, we insert a Bose-Einstein distribution 
function. For comparison we also consider the case of a 
Maxwell-Boltzmann distribution. Then we obtain 



ei(IPl) 



c vac 
"1 



2T 

l + -p^X 

HP 



/ 1 - exp (- 
In V 



2T 



1 — cxp ^~ 



2T 



-2ft 



BE , 

MB . 

(32) 

The resulting expression depends on time t = Xq via tem- 
perature T and chemical potential fi of the toy baryons. 
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For the rest of this section we assume <C T, as in real- 
istic scenarios of leptogenesis, for the purpose of illustra- 
tion 11 . The temperature and momentum dependence of 
the medium correction in the range of typical momenta 
|p| ~ T is shown in the shaded areas in Fig. [5] for the 
BE and MB cases, respectively. It is instructive to con- 
sider the nonrelativistic (NR) regime (T, |p| <C Mi) and 
the ultrarelativistic (UR) regime (T > |p| ^> Mi). In the 
nonrelativistic (NR) limit, the BE and MB cases coincide, 
and the medium correction is exponentially suppressed, 



ei(|p|)/er ^1 + 2 exp 



MA 
2T J 



(33) 



Furthermore, it is independent of the momentum |p|. In 
the UR limit, the medium correction for the BE and 
MB cases behaves quite differently: In the MB case the 
medium correction saturates at ei/e'[ ac < 3. In the BE 
case, it is logarithmically enhanced 12 (sec Fig. [5]), 



ei(|p|)A? e 



IT 
r|p| 



In 



4T|p| 



Ml 



8p 2 m 2 
Afj'Cl+r) 



(34) 



This effect is due to Bose enhancement. Thus, we find 
that the quantum statistics is important for the medium 
correction. In the following section, we will see that the 
logarithmic enhancement at high energies is also sup- 
pressed by the inclusion of sizable negative chemical po- 
tentials (which is necessary within the toy model, see be- 
low). In Fig.0 we also show the CP-violating parameter 
(ei) obtained from averaging Eq. (j32[) over the momen- 
tum |p| (see AppendixlFjl . As expected, (ei) ~ ei(|p| ~ 
T). 

Before discussing the impact of the medium correc- 
tion quantitatively, we would like to comment on the re- 
lation between the Kadanoff-Baym (top-down) and the 
canonical (bottom-up) approach. As has been men- 
tioned before, in vacuum, the top-down result Eq. ([25]) 
for the CP- violating parameter coincides with the canon- 
ical result Eq. (j2j). Nevertheless, we emphasize again 
that the structure of the Boltzmann equations differs be- 
tween the two approaches, i.e. the former are free of the 
double-counting problem. Furthermore, it is also impor- 
tant to note that the size of the medium correction dif- 
fers between the top-down and the bottom-up approach. 



11 In realistic scenarios, the leptons and Higgs fields arc in equi- 
librium with gauge bosons such that fi = —p,. The smallness of 
the asymmetry then ensures the smallness of the chemical poten- 
tials, |/i|/T = \p\/T ~ 10~ 10 . Within the toy model, due to the 
absence of gauge interactions, it is possible to have \fi\ ~ \p\ ~ T 
while the asymmetry remains small. It turns out that this is 
even necessary to obtain consistent numerical solutions within 
the present scenario, see Sec. IIVI 

12 Note that the logarithmic enhancement at high energies is cut-off 
for extremely high energies (|p| S> M^/m) in the UV due to the 
second summand of the denominator inside the logarithm. Since 
we assume that m <C Afi/10, we can neglect this term in the 
relevant temperature range M\/T > 0.1. 



Within the latter, the medium corrections have been dis- 
cussed by replacing e" ac — > e| in the canonical Boltz- 
mann equations. Hereby ef L involves the vertex loop cal- 
culated within thermal field theory (see e.g. fl9l.l20l]). For 
the toy model, e* is given in Eq. (|A5[) of AppendixlAl It 
involves an additional term compared to the top-down 
result (|28[). which is quadratic in the particle distribu- 
tion function. In Fig.|5]we show that the medium correc- 
tion would be significantly over-estimated in the canoni- 
cal thermal field theory approach within the toy model. 

For realistic models of leptogenesis, the vertex loop 
contains scalar and fermionic lines in general. In con- 
trast to the scalars, the latter tend to decrease the size 
of the CP-violating parameter. Therefore the results 
shown here can only be used indirectly to make state- 
ments for phenomenology. For the standard scenario of 
thermal leptogenesis with hierarchical right-handed neu- 
trino spectrum, it has been observed in [l9( that, within 
thermal field theory, the effects of Bose enhancement and 
Pauli blocking tend to cancel each other. However, since 
the medium corrections differ within the Kadanoff-Baym 
formalism, this cancellation may no longer occur. There- 
fore one might expect that the medium correction is un- 
derestimated in this case, contrary to the situation en- 
countered within the toy model discussed here. 

Note that we have neglected the effect of thermal 
masses here for simplicity. Their impact on leptogenesis 
has been studied e.g. in [2fj] within the framework of ther- 
mal field theory. For the toy model considered here, it is 
consistent to neglect thermal masses due to the absence 
of gauge interactions. We stress, however, that it is also 
possible to include time-dependent effective masses sys- 
tematically within the Kadanoff-Baym formalism. One 
of the effects caused by the thermal masses of the com- 
plex field is the cutoff of the logarithmic enhancement 
of the CP- violating parameter at high energies. In ad- 
dition, in an asymmetric medium the effective masses of 
the particles and antiparticles are not equal, which leads 
to an effective CP violation. This effect is studied in 
[HBI . The medium corrections to the masses of the heavy 
real fields can be important in the case of the degenerate 
mass spectrum, which we investigate in [29j . 

Since the quantum statistic is important for the CP- 
violating parameter in medium, we also have to include 
the quantum statistical terms in the gain and loss terms 
of the Boltzmann equations for consistency. This will be 
discussed in the following section. 



IV. NUMERICAL RESULTS 

To obtain the three Boltzmann equations for /, / and 
f^ ± , we integrate each of Eqs. (fl4|) . (fl~5|) and (|20|) over 
Po (left- and right-hand side) and choose the positive en- 
ergy solution [30(. In agreement with the cosmological 
principle, we solve the system of Boltzmann equations 
with spatially homogeneous and momentum isotropic 
distribution functions in (flat and radiation dominated) 
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Friedman-Robertson-Walker (FRW) space-time. In this 
case the left-hand side of the Boltzmann equation is given 
by 

L[/](|p|) =p a X> a /(|p|) 

= p0 {i-^ H W\) mh (35) 

where H = a/a is the Hubble parameter. Since in 
the quasiparticlc approximation the spectral functions 



are proportional to 8(p 2 — m 2 ), the integration over the 
time components of the quasiparticlc's four-momenta can 
be performed trivially in each of the corresponding self- 
energies. After the integration the volume element dH p 
is replaced by dll 3 = d 3 p/(27r) 3 /2E, where E is energy 
of the on-shell quasiparticle. 



With these modifications, the network of quantum- 
corrected Boltzmann equations takes the form: 



L[f}(\Pl) = CW4/KIP|) + C bb ^ bb [f, /](|p|) + CW^J/, jVJGpl) , (36a) 

L[f](\p\) = c- bb „ bb [f}(\p\) + Cumuli, M\p\) + Cu^ALUAM) , (36b) 

^ I ](|p|) = C^^[^ 1 ,/](|p|)+C v , 1 ^[/^ 1 ,/](|p|), (36c) 

where the different collision terms for the transition between two states i and / are denoted by Ci++f. Note that, due 
to the isotropy of the FRW universe, the collision terms also depend only on the absolute value of the momenta. For 
the 2 — 2 scattering processes in (|36a|) we find 

c bb ^ bb [f](\ P \) = hJdnl 2 dnl 3 dnl 4 (27T) 4 s^(p + P2-P3-P4) 

x iA 2 {[l + /(|p|)][l + /(| P2 |)]/(|P3|)/(|P4|) - /(|p|)/(| P2 |)[l + /(| P3 |)][1 + /(|P4|)]} , (37a) 

^ 65 [/,/](|p|) = Udnl 2 dnl 3 dnl 4 (27T)^(p + P2 - P3 - Pi ) 

x A 2 {[1 + /(|p|)][l + /(|P2|)]/(|P3|)/(|P4|) - /(|P|)/(| P2 |)[1 + /(| P3 |)][1 + /(|P4|)]} • (37b) 

The corresponding terms in the equation for 6 can be obtained by replacing / with / in (|37a[) and (|37b[) . If the 
generated asymmetry is small, as we assume here, then / ~ /. In this case the CP-violating contributions to the 
right-hand side of (|36cp cancel out and we obtain: 



^^[/^JKIpD + C^^I^JKIpI) ~ \^dlil 2 dlil 3 {2^\p-p 2 -p 3 ) 

xi|.9i| 2 ({[l + /^(|p|)]/(| P2 |)/(|P3|) - /v,dPl)[l + /(|P 2 |)][1 + /(IPsI)]} 

+{[1 + M(|P|)]/(|P2|)/(|P3|) - /*(|p|)[l + /(|P 2 |)][1 + /(|P3|)]}) ■ (38) 

The collision terms for the (inverse) decay of the heavy particle into a bb or bb pair explicitly contain the CP- violating 
parameter e defined in Eq. (|31[) : 

cw 01 \f, /^](| P |) = i / dul 2 dnl 3 {2*Y8^(p 2 -p- P3 ) 

x |.9i| 2 [l + ex(|p 2 |)]{[l + /(|p|)][l + /(| P3 |)]/ 01 (| P2 |) - /(|p|)/(|p 3 |)[l + JW(|p 2 |)]} (39a) 
Cu^ [/, iW](|p|) = \ J dU p2 dU p3 (27:)^(p 2 -p-pa) 

x |.9i| 2 [l - ei (|p 2 |)]{[l + /(|p|)][l + /(|p 3 |)]/^(|p 3 |) - /(|P|)/(|P 3 |)[1 + UM)}} ■ (39b) 



The factors 1/2 associated with the couplings in ()37ap 
and (|38|) correctly account for the symmetrization of col- 
lision integrals which include integration over the mo- 
menta of two identical particles in the initial/final state. 
They have been consistently obtained in the derivation 
from the Kadanoff-Baym equations. We would like to 
stress again that the structure of (f3T)j) differs from the 
usual structure obtained in the conventional bottom-up 
approach. In particular, we do not need to include col- 



lision terms for the processes bb -O- bb (not even the RIS 
part of it) because our collision terms for the processes 
bb O ipi and bb •H- tpi do not suffer from the generation of 
an asymmetry in equilibrium. The network of Boltzmann 
equations (|36[) should be understood in the generalized 
sense: the "amplitudes" differ from the usual perturba- 
tive matrix elements and do not have their symmetry 
properties. 

To stay consistent in our model we also keep the quan- 
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turn statistical factors for bosons. The implications of 
this new structure for a phcnomcnological theory of lcp- 
togcncsis will be discussed elsewhere. To study the ef- 
fect of the quantum corrections, we compare the results 
obtained by integrating the network of Boltzmann equa- 
tions with quantum-corrected ei(|p|) to those which are 
obtained after replacing ei(|p|) with e\ ac . This means 
that we keep here the new structure of the Boltzmann 
equations and study corrections which arise from the 
quantum-corrected t\ only. In the vacuum limit the 
structure of Eq. ([36)) corresponds to the one which has 
been assumed implicitly in [56| if the quantum statisti- 
cal factors and the symmetrization factor are replaced 
accordingly. 

In both cases we start at sufficiently high tempera- 
tures so that all species, including with mass M\ = 
10 10 GeV, have relativistic initial abundances which cor- 
responds to the most frequently discussed case. 13 Be- 
cause of the presence of the statistical factors we need to 
start with sufficiently negative chemical potentials as to 
avoid Bose-Einstein condensation of the different species 
during their evolution. 14 We choose them such that they 
are related by = 2/i^ = 2/xg, i.e. the system is in 
chemical equilibrium. 

The coupling A can be adjusted such that the rates of 
the 2 — 2 interactions ([57)) are much larger than those 
of the decays and inverse decays (|38| - (I5^1) at all times 15 . 
This keeps b and b close to kinetic equilibrium, just as 
Higgs particles and leptons are kept in equilibrium by 
rapid gauge interactions in the standard scenario. The 
distribution functions for these species are therefore given 
by their equilibrium form throughout the entire evolu- 
tion. This means that they can be described in terms 
of four parameters /j,b, Tb and /i^, T5. The interactions 
(J37J enforce the relation Tj, = Tb between the parame- 
ters. Therefore, it is sufficient to study the evolution of 
/ and / in terms of the remaining three parameters (see 
Appcndix lGj) . The evolution of ipx, however, is studied 
in terms of the complete distribution function discretized 
on a grid with 400 momentum modes. Our computation, 
therefore, includes classical nonequilibrium effects in the 



2:2 



13 Another scenario, which is frequently discussed in the literature, 
is that the Majorana neutrinos could have zero initial abundance. 
In this case we would expect the differences in the time evolu- 
tion of the asymmetry to be larger in general. However, this can 
have an effect on the final asymmetry only if the asymmetry, pro- 
duced in the intermediate step of thcrmalization, is not washed 
out again before the Majorana neutrinos decouple, i.e. for small 
washout factor K < 1. 

14 In this regime it would not be appropriate to describe the sys- 
tem by conventional Boltzmann kinetic equations. Since we are 
interested in scenarios that are qualitatively similar to realistic 
models of leptogenesis here, we do not consider this case. 

15 As we will show in AppcndixlGl in this case there is no need to 
compute the collision integrals for 2 — 2 scattering explicitly and 
we can use perturbative values for A for most of the relevant 
range of \g\ 2 assuming that it is sufficient to demand that the 
rate for bb <-> bb is at least 10 3 times larger than that of bb -f> ipi. 
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FIG. 6: The ratio (ei) /ef 1 ". The shape of the curves differs 
from that of the corresponding graph in Fig. [5] mainly because 
its computation involves a finite chemical potential (which 
depends on Mi/T) here. Similar graphs can be obtained by 
including a finite chemical potential in Eq. (I32[) . 



decay of ipi. Such effects have been studied recently in 
[5a - l58l ] . All integrals are evaluated numerically including 
all quantum statistical factors for stimulated emission. 
We define the generated "baryon" asymmetry as 



r)(Mi/T) 



s(M x /T) 



(40) 



Here rib and n^, the number densities of species b and b 
respectively (compare (|G7[) h are computed in the pres- 
ence of the quantum-corrected ei, and s is the standard 
entropy density [l| . The analogous asymmetry computed 
with e\ ac is denoted by rf ac {Mi/T), 

Figure [6] shows the numerical value of the ratio 
(ei) /e\ ac for various values of the washout parameter 
K = T/H{M 1 ) = |.gi| 2 m pi /(4.5-167r v /^Af 1 3 ). The flatten- 
ing for small M\/T as compared to the thermal equilib- 
rium result in Fig.[5]is due to the finite chemical potential 
of b. This shows that larger corrections could be obtained 
if additional interactions for b and b are introduced which 
would allow one to start with smaller chemical potentials 
and hence lead to a stronger enhancement. 

The buildup of the asymmetry with and without quan- 
tum corrections as a function of the inverse temperature 
is depicted in Figs. [7] and [5J Comparing these figures one 
can verify the enhancement of the asymmetry at interme- 
diate stages for larger washout factors (case d) . Note also 
that due to the medium contribution to the CP-violating 
parameter the generated asymmetry is not a monotonous 
function of the washout parameter k. 

The dependence of the resulting final asymmetries 
r) = r/(Mi/T -> 00) and rj vac = j ] vac {M 1 lT -> 00) as well 
as the ratio (77 — r j vac y r j vac on the washout parameter is 
presented in Fig.[SJ The asymmetry is always larger when 
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FIG. 7: The asymmetry r)(Mi/T) with quantum corrections 
included. In the weak washout regime (case a) the asymmetry 
is produced at smaller temperatures and it is not necessarily 
larger than for larger washout factors (compare a and b). 
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FIG. 9: The final asymmetries and the ratio (rj - r) vac )/r) vac 
over washout factor k. The cases a, b, c, d, e, f correspond 
to washout factors 0.01, 0.1, 0.366, 1, 10, 100. Case c is close 
to the maximum relative excess of the quantum-corrected re- 
sults at k — 0.34. In contrast to the usual results the final 
asymmetry does not take its maximum value for the smallest 
washout factor. Instead, the asymmetry 77 peaks at k ~ 0.059. 



creasingly dominates E\ = (M-j 2 + p 2 )^ as the momenta 
are redshifted to smaller values. This means that the 
integration is over an interval in which the distribution 
/ becomes smaller and smaller. This explains why the 
quantum corrections tend to zero for small k. For the 
same reasons the relative effect of the quantum correc- 
tions peaks at a moderate k ~ 0.34 with about 26%. 

We note again that the size and the sign of the correc- 
tions depend on the quantum statistics of the particles 
in the vertex loop and will be different in a phenomcno- 
logical scenario. Further plots and details about the nu- 
merical algorithm can be found in AppcndixlGl 



FIG. 8: The asymmetry rf ac (M^lT) without quantum cor- 
rections. 

quantum corrections are taken into account compared to 
the results without corrections (compare Sec. lIII|) . The 
asymmetry 77 has a maximum for moderate washout fac- 
tors k ~ 0.059 in contrast to the usual result which has 
its maximum in the limit of zero washout factor. Our in- 
terpretation of this result is as follows: For large washout 
factors the enhancement of e\ due to the quantum correc- 
tions enhances the asymmetry generated by the decays 
only at intermediate stages, because the same processes 
diminish the asymmetry in particular at late times where 
the averaged asymmetry drops to smaller values (com- 
pare Fig.|B]). For small k the particles decay late, and the 
backreaction is largely suppressed so that the washout is 
ineffective. However the interval of integration in Eq. (|31j) 
is located at relatively large momenta since the mass in- 



V. CONCLUSIONS AND OUTLOOK 

In this paper, we have studied leptogenesis in a 
simple toy model consisting of one complex and two 
real scalar fields in a top-down approach, using the 
Schwinger-Keldysh/Kadanoff-Baym formalism as the 
starting point. This treatment, based on noncquilib- 
rium quantum field theory techniques, is motivated by 
the fact that it allows a unified description of two key in- 
gredients of leptogenesis, namely deviation from thermal 
equilibrium and loop-induced CP violation. It has sev- 
eral important advantages in comparison to the canoni- 
cal bottom-up (Boltzmann) approach. In particular, the 
full Kadanoff-Baym equations do not rely on the con- 
cept of quasiparticles and their collisions in the plasma. 
However, if the quasiparticle picture is applicable, as we 
have assumed here, the Kadanoff-Baym formalism con- 
sistently accounts for the dependence of the quasiparti- 
cles' properties as well as scattering and decay rates on 
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the state of the medium. 

The out-of-equilibrium dynamics of the quasiparticles 
is described by a system of approximate self-consistent 
kinetic equations - quantum- corrected Boltzmann equa- 
tions - which we have derived here starting from the full 
system of the Kadanoff-Baym equations. We find that 
the structure of the quantum-corrected Boltzmann equa- 
tions automatically ensures that no asymmetry is pro- 
duced in thermal equilibrium. In other words there is no 
need for the real intermediate state subtraction, i.e. the 
formalism is free of the double-counting problem typical 
for the canonical approach. 

One of the key quantities in leptogenesis is the CP- 
violating parameter. Earlier studies have shown that 
there are two sources of CP violation: self-energy and 
vertex contributions. In this work, we have concentrated 
on the latter one. We have found that for scalar fields 
the medium effects increase the vertex contribution to 
the CP-violating parameter. At high temperatures it is 
up to an order of magnitude larger than in vacuum and 
asymptotically approaches the vacuum value as the tem- 
perature decreases. This result can be traced back to a 
Bose enhancement of the vertex loop correction. In the 
Maxwell-Boltzmann approximation, the corresponding 
CP-violating parameter is increased at most by a factor 
two. We would also like to note that, in the vacuum limit, 
the CP-violating parameter obtained via the Kadanoff- 
Baym formalism agrees with the value obtained within 
the canonical formalism, as expected. 

It is interesting that, contrary to the results obtained 
earlier in the framework of thermal field theory by replac- 
ing the zero-temperature propagators with finite temper- 
ature propagators in the matrix elements of the Boltz- 
mann equation, the medium corrections depend only lin- 
early on the particle number densities. Stated differently, 
only one of the internal lines in the vertex loop is "ther- 
mal". Moreover, the medium corrections to the vertex 
CP-violating parameter depend only on the density of 
the toy baryons and are independent of the density of 
the "Majorana" particles. Since the decaying heavy par- 
ticles as well as the intermediate on-shcll states propagate 
with respect to the thermal bath's rest frame, the CP- 
violating parameter in each individual decay depends on 
the momenta of the initial and final states. 

We have solved the system of the quantum- corrected 
Boltzmann equations numerically. Because of the 
medium corrections to the CP-violating parameter the 
asymmetry reaches its maximum value at a small but fi- 
nite value of the washout parameter n, rather than for 
k — v 0, as it is the case in the canonical approach. To 
avoid the regime of Bose-Einstein condensation we have 
to assume that the species have rather large chemical po- 
tentials initially. This decreases the medium correction 
to the CP-violating parameter. As a result the generated 
asymmetry differs from its value in the canonical formal- 
ism by approximately 26%. However, in a scenario in 
which the chemical potentials are close to zero the devia- 
tion could reach the 100% level. As has been mentioned 



above, our results differ from the results of the calcula- 
tions performed in vacuum and in the framework of ther- 
mal field theory. Therefore, we argue that one should use 
the quantum-corrected Boltzmann equations (or the full 
Kadanoff-Baym equations). On the other hand, to ob- 
tain order of magnitude approximations, it seems safe to 
use the canonical approach in the present scenario. 

The techniques described in this work can also be 
used to study quantum noncquilibrium effects within 
phcnomcnological scenarios of leptogenesis. In particu- 
lar, the technical advantages of the Kadanoff-Baym for- 
malism demonstrated above, like the absence of double- 
counting problems, are quite generic and therefore should 
not depend on the details of the model. Furthermore, we 
expect that the difference in the size of the medium cor- 
rections, obtained from the Kadanoff-Baym formalism 
and thermal field theory respectively, persists if our ap- 
proach is applied to such scenarios. 

We would like to stress again that, in this paper, we 
have considered only the vertex contribution to the CP- 
violating parameter. The self-energy contribution is com- 
parable to the vertex contribution or, in the case of res- 
onant leptogenesis, considerably larger than the vertex 
contribution. We will study it in a forthcoming pa- 
per [2^]. In addition in [55[ we plan to investigate the 
influence of a nonzero asymmetry on the CP-violating 
parameter. Furthermore, for phcnomcnological models, 
thermal masses can become relevant. These can also 
consistently be described within the Kadanoff-Baym for- 
malism. Finally, we note that it would also be inter- 
esting to investigate numerical solutions of the full set 
of Kadanoff-Baym equations without further approxima- 
tions as some of their properties cannot be included in 
Boltzmann- like equations for principle reasons. This is 
however beyond the scope of the present work. 
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Appendix A: CP-violating parameter in the 
bottom-up approach 

In this appendix, we review the calculation of the ver- 
tex contribution to the CP-violating parameter in vac- 
uum, e\ ac = - T^.^n) I (Tfa-tbb + r^^bft), in 
the conventional in-out formalism. It is generated by the 
interference of the tree-level and one-loop amplitudes (see 
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Fig.®, 



M 



(0) _ _ . * 
ipi^bb ~ l 9i 

<(1) 



M^Ubb = ^9*9* m?C (M?, 0, Mj) . 

In the limit of massless toy baryons, the scalar 1-loop 
vertex three-point function Cq is given by [EH |60| : 



C Q ( P lplMf) = 



(iTr 2 )- 1 ^ 



(q 2 - Mf)(q +p 2 ) 2 (q + p 2 ~ PlY 



(A2) 



where p\ and p 2 are momenta of the decaying heavy par- 
ticle and of one of the decay products, respectively. The 
tree-level and one-loop amplitudes of the decay process 
%jj — > bb differ from (|A1|) only by conjugation of the cou- 
plings. Therefore, at leading order, we obtain for the 
CP-violating parameter: 



Im 



9i9i 



ImC (M?,0,M?). (A3) 



8tt 2 \9l9j. 
Substituting the result for the three-point function, 



C o (Mf,0,M, 2 



1 



Li 2 




6 



. (A4) 



(with dilogarithm Li2 as defined in [60j]) into Eq. (|A3[) . we 
obtain Eq. @ . The same calculation can be performed 
within thermal quantum field theory by using thermal 
propagators in Eq. (|A2[) . The result is [l9l. [53l| 



8tt M? 



-Im 



5,: 5j 



dn 1 + /f + /* 



th fth 



2/r/ : 



4tt M?/M? + |(1 + cos( 



(A5) 



where f} h = [exp((E ki - /x)/T) - l]" 1 , see Eq. ([29]). The 
ellipsis denote similar contributions involving the distri- 
bution function of the heavy scalar. These are suppressed 
in the hierarchical limit. 



Appendix B: Kadanoff— Baym formalism for the 
complex scalar field 

Here, we derive the Kadanoff-Baym and quantum- 
corrected Boltzmann equations for the complex scalar 
field. 



1. Schwinger— Dyson equation 

Our starting point is the generating functional for 
Green's functions [oil ]: 

Z[J, K] = J @b3>bexp[i{S + Jb + Jb + bKb)} , (Bl) 



FIG. 10: Closed real-time path C. 



where the field and the external sources are defined on 
the the positive and negative branches of the Schwinger- 
Keldysh closed real-time contour shown in Fig. [T(Ill3ll .l62l- 
lo6| . The scalar products of the local and bilocal sources 
J(x) and K(x, y) and the field are defined as invariant 
configuration space integrals [30. l67ll . Furthermore, we 
use the compact notation of Ref. 31[ for contour integrals 
over the closed real-time path. Note that the sources are 
now complex functions. The requirement that the last 
term in (JbTJ) be real implies that K(x,y) = K*(y,x). 

The functional derivatives of the generating functional 
for connected Green's functions, 



W[J,K] = -ilnZ[J,K] , 
with respect to the external sources read 



dW[J,K] 

dJ(x) 
dW[J,K] 
dK(x, y) 



B{x) 



= \[D{y,x) + B{x)B{y)] 



(B2) 

(B3a) 
(B3b) 



B and D denote the expectation value and the propaga- 
tor of the field respectively. The derivative of W with 
respect to J is just the complex conjugate of (|B3aj) . 

Performing a Legendre transform of the generating 
functional for connected Green's functions, we obtain the 
effective action 

T[D,B}= W{J, K]- JB-JB- tr[KD] - BKB. (B4) 



Making use of the chain rule and Eqs. (|B3|) . we find for 
functional derivatives of the effective action 



6T[D,B] 

5B(x) 
ST[D, B] 
6D(x,y) 



= -J(x) - f@ 4 z K(x,z)B{z) , (B5a) 
= -K(y,x). (B5b) 



Next, we shift the complex field by its expectation 
value b — > b + B. Since the integration measure in the 
path integral is translationally invariant, the effective ac- 
tion can be rewritten in the form 

T[D, B] = - i In J %>b3ibciip[i{S + Jb + Jb + bKb)} 

+ Sd[B]-ix[KD]. (B6) 

Now, we tentatively write the effective action in the 
form [6l| 

T[D,B] = S c i[B] -M bidet [D^ 1 ] +itr \_9~ X D\ 

+ T 2 [D,B], (B7) 
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thus defining the functional T 2 . 

The third term on the right-hand side of (|B7[) is given 
by a convolution of the field propagator D and the free 
inverse propagator Its differentiation with respect 

to D(y,x) gives [30l ] 



- 1 (x,y)=i(n x +m 2 )6°(x,y) 



(B8) 



The functional derivative of the second term on the right- 
hand side of (|B7j) can be obtained upon use of 

f2> i zD- 1 (x,z)D(z,y) = 53(x 1 y), (B9) 

and is given by — iD~ l {y, x) . Consequently, we obtain 

5T\D,B] . „ . 5T 2 \D,B} 

1 ' J = -iD~ 1 (y,x) +i@- 1 (y,x) + 1 ' J 



6D(x,y) 



K(y,x) . 



5D(x,y) 

(BIO) 



Physical reality corresponds to vanishing sources. 16 
Therefore, Eq. (|B10[) can be rewritten in the form 



D- 1 (x,y) = ®- 1 (x,y)-Ti(x,y). 
where the self-energy is defined by 

.ST 2 [D,B] 



T,(x,y) 



SD(y,x) 



(Bll) 



(B12) 



Note that the factor two in the definition of the self- 
energy [3(| E3] is absent, just as one would expect for a 
complex field. 



2. Kadanoff— Baym equations 



Convolving the Schwinger-Dyson equation (jBll[) with 
D from the right and using Eq. (|B9[) . we obtain 



i[D x + m 2 ]D(x, y) = 6 g (x, y) 

+ / @ 4 zT,(x,z)D(z,y) . (B13) 

Next, following the usual procedure, we represent the 
time-ordered propagator as a linear combination of the 
statistical propagator and spectral function: 

D(x,y) = D F (x,y) - |sign c (x° - y°)D p (x,y) , (B14) 

where sign c denotes the signum function with respect to 
time ordering along the closed time path, and 

D F (x,y) = \{[b(x), %)]+), (B15a) 
D p (x,y) = i([b(x),b(y)}-), (B15b) 



16 To be precise, within nonequilibrium field theory, this is only true 
for times x°,y° > The local and bi-local sources supported 

at x° = y° = ti n it formally encode the information about the 
(Gaussian) initial state (see e.g. [HI])- However, these sources 
do not appear explicitly in the Kadanoff-Baym equations, and 
therefore we omit them here. 



where the subscripts "+" and "— " denote the anticom- 
mutator and the commutator of the fields. 

To find out how Dp and D p behave under complex 
conjugation let us introduce 

£>> (x, y) = (b(x)b(y)) = Tr[^ b(x)b(y)} , (B16a) 
D < (x, y) = (b(y)b(x)) = Tr[0> b(y)b{x)] . (B16b) 

Using the Hcrmiticity of the density matrix and the 
cyclic invariance of the trace, we obtain 

D* > (x,y) = D > (y,x), £>* (x, y) = £>< {y, x) . (B17) 

Consequently 

D F (x,y)=D F (y,x), D* p (x, y) = -D p (y, x) . (B18) 

Analogous relations also hold for the spectral and statis- 
tical components of the self-energy. 

The local part of the self-energy is proportional to the 
Dirac 5 function and can be absorbed in the effective 
mass of the field, m 2 (x) = m 2 + Y} oc (x,x) , whereas the 
remaining part of the self-energy can be split into a spec- 
tral part Yi p and a statistical part Tip in a complete anal- 
ogy to (|B14l) . 

Because of the sign function, the action of the Laplace- 
Beltrami operator on (|B14|) gives rise to a product of 
g 00 5(x° ,y°) and Vgf p (a;,y). Upon use of the definition 
of the spectral function and canonical commutation re- 
lations for a complex scalar field this product reduces to 
the generalized S function S 9 (x, y) j3(|, which cancels the 
5 function on the right-hand side of (|B13|) . 

Separating spectral and statistical components in 
Eq. (|B13[) . we obtain a system of Kadanoff-Baym equa- 
tions very similar to that for the real scalar field (30l |: 

y° 

[D x + m 2 {x)]D F {x, y) = J @ 4 z T,p(x, z)D p {z, y) 

o 

x° 

- f@ i zTi p (x,z)D F (z,y) : (B19) 
o 

y° 

[a x +m 2 (x)]D p (x,y) - J^z E p (x, z)D p (z, y) . (B20) 

x° 

One should, however, keep in mind that the functions 
in (|B19[) and (|B20[) are complex. That is, we get four 
equations for the real and imaginary components of the 
spectral function and the statistical propagator. The 
complete information about the (Gaussian) initial state 
specified at the initial time £,„,j = enters via the ini- 
tial conditions of the two-point functions Dp, d x oDp, 
dyoDp and d x od y aDp evaluated at x° — y° — t init . The 
corresponding initial conditions for the spectral function 
are fixed by the equal-time commutation relation of the 
complex field, as for the real case (38l |. 

We note that the proper renormalization of the 
Kadanoff-Baym equations (|B19IB20p would require one 
to also take non-Gaussian correlations of the initial 
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state into account [68|, |69(. However, the derivation of 
quantum-corrected Boltzmann equations considered here 
involves the limit t,- n ,4 — > — oo, and therefore the initial 
correlations should have a negligible effect. 

3. Quantum kinetics 

The Kadanoff-Baym equation for the statistical prop- 
agator (spectral function) can be rewritten in terms of 
the advanced and retarded propagators, D R and D A : 

[U x + m 2 {x)]D F{p) {x,y) = - J @*z9(z°) 
x [Z F{p) (x,z)D A (z,y) + Z R (x,z)D F(p) {z,y)}. (B21) 

Because of (|B18[) . the retarded and advanced propagators 
are related by 

D R (x,y) = 6(x°-y°)D p (x,y) 

= -0(x° - y°)D* p (y, x) = D* A (y, x) . (B22) 

Therefore, after interchange of x and y in (|B21[) and com- 
plex conjugation of the resulting equation, we find 

[n y +m 2 (y)]D F{p) (x,y) = -J^z0(z°) 
x [Dr(x, z)Y, F{p) (z,y) + D F{p) (x, z)Y, A (z,y)} . (B23) 

Since (|B23|) has been obtained from (|B21|) by reversible 
transformations, a solution of (|B21[) is also a solution 
of (|B23I) . Consequently, it is also a solution of the sum 
(which will be referred to as constraint equation) and the 
difference (which will be referred to as kinetic equation) 
of (1B2TT) and £B23j). 

To analyze the constraint and kinetic equations it is 
convenient to introduce the center and relative coordi- 
nates, X and s [5l[. In terms of the center and relative 
coordinates relations (|B18j) can be rewritten in the form 

D* F (X, s) = D F (X, -s), D* p (X, s) = -D P (X, -s) . 

Consequently, even in the case of a complex scalar field, 
the Wigner transforms of the spectral function and sta- 
tistical propagator, 

D F (X,p) = ^/ = <? x J d 4 s e ips D F (X, s) , (B24a) 
D p (X,p) = -iV=9xJ d A se^D p {X,s) , (B24b) 

are real-valued functions. The Wigner transforms of the 
retarded and advanced propagators are defined analo- 
gously to (|B24ap . From (|B22I) it then follows that rela- 
tion 

D A {X,p)=D R {X,p) (B25) 

also holds for a complex scalar field. Another very useful 
relation, 

D R (X,p)-D A (X,p)=W p (X,p), (B26) 



results from the definitions (|B22[) and (|B24|) and the 
equality 8(s°)+6(~s°) = 1. Equations (jB25)) and (|B26| 
in particular imply that 

D R(A) (X,p) = D h (X,p) ± lD p {X,p), (B27) 

where Dh(X,p) = ReD R (X,p) has been introduced. A 
similar relation also holds for the retarded and advanced 
self- energies. 

Let us now subtract (|B23[) from (|B21j) and then Wigner 
transform the left- and right-hand sides of the resulting 
equation. Furthermore, we send the initial time to the 
infinite past, ti n n — > — oo, i.e. we drop the functions 9(z°) 
on the right-hand sides of Eqs. ()B21|) and (|B23|) . Phys- 
ically, this means that we neglect the effects of initial 
correlations. Additionally, we perform a gradient expan- 
sion with respect to the central coordinate X. Proceeding 
as in [3^1, we obtain a kinetic equation for the spectral 
function. To linear order in the gradients it reads 

{uj(X,p),D p (X,p)} PB = {X p (X,p),D h (X,p)} PB , (B28) 

where we have introduced 

uj(X, P ) = g^p^pu - m 2 (X) - V h {X,p) , (B29) 

and the Poisson brackets are defined by [3(3] 

{A(X,p),B(X,p)} PB = ^A(X,p)V a B(X,p) 

-V a A(X,p)^-B(X,p). (B30) 
op a 

Wigner transforming the sum of (|B23|) and ()B2ip . we 
obtain the constraint equation for the spectral function. 
To linear order in the gradients it is an algebraic equation: 

u(X,p)D p (X,p) = X p (X,p)D h (X,p) . (B31) 

To close the system and to analyze the spectrum, we 
also need the equations for the retarded and advanced 
propagators. They can be obtained from (|B21[) and (|E>23|) 
upon use of the definitions of D R and D A and the canon- 
ical commutation relations: 

[D x + m 2 (x)]D R{A) (x, y) = S 9 (x, y) 

- J ® A z E fl(A) (x, z)D R{A) (z, y) , (B32) 
[□„ + m 2 (y)]D A(R) (x, y) = S 9 (x, y) 

-f@ 4 z D A(R) (x, z)Z A(R) (z,y). (B33) 

Wigner transforming the difference of (|B32j) and (IB33|) 
and subtracting (jB28|) . we obtain the kinetic equation for 
real part of the retarded and advanced propagators: 

{u(X,p),D h (X,p)}p B = 

-i{S p (A,p),^ p (A,p)} PB . (B34) 

Wigner transforming the sum of (|B32|1 and (|B33[) and 
subtracting (|B31|) . we obtain the second constraint equa- 
tion: 

uj{X,p)D h {X,p) = -1 - \V p {X,p)D p {X,p) . (B35) 
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The solution of the system of constraint equations (IB31|) 
and (|B35|) reads 



D p (X,p) 
D h (X,p) 



u(X,p) 



S P (X,p) 



D p (X,p). 



(B36a) 
(B36b) 



As can be checked by substitution, solution (|B36[) is also 
solution of the kinetic equations (|B28j) and ()B34jl . In 
other words, to linear order in the gradients we have an- 
alytic expressions for the spectral function and retarded 
(advanced) propagators. The spectral function has a 
sharp peak on the mass shell, i.e. for us(X,p) = 0. The 
height and exact shape of the peak are time-dependent. 
Proceeding in a similar way, we can derive the kinetic 



Z F (X,p)D p (X,p) , (B37) 



{o J (X,p),D F (X,p)} PB - 
+ D F (X,p)Z p (X,p) 

and the constraint, 

u(X,p)D F (X,p) = \{Z F (X,p),D p (X,p)} PB (B38) 
+ \{D F (X,p), Z p (X,p)} PB + E F (X,p)D h (X,p) 

equations for the statistical propagator. The constraint 
equation is no longer algebraic and cannot be, generally 
speaking, solved analytically. However, if the system is 
in thermal equilibrium, then all the quantities are con- 
stant in time and space and the Poisson brackets in (IB38|) 
vanish identically. The solution of the resulting equation 
reads 



(B39) 



That is, we have obtained the fluctuation-dissipation 

As can 



relation from the constraint equation (IB38|) 
be checked by substitution, in equilibrium (|B39j) is in- 
deed a solution of (|B37[) . Furthermore, using (|B14[) and 
the Kubo-Martin-Schwinger periodicity condition we find 

D e F q ( P )= [f eq (p)+ \]D e p q {p), (B40) 
where f^ eqS> is the Bose-Einstein distribution function. 

4. Quantum-corrected Boltzmann equations 

The spectral function for the complex scalar field 
(|B36a|) has a Breit-Wi gner shape and peaks at lo = 0. 
The height of the peak is inversely proportional to the 
spectral self-energy and tends to infinity in the limit of 
vanishing coupling constant. Since, furthermore, the area 
under D p (X,p) is constant [30| to a first approximation 
it can be replaced by a Dirac 8 function: 



D p (X,p) = 27rsign(p ) 6 {g pv p»p v 



m 



(B41) 



Equation (|B41|) is referred to as quasiparticle approxima- 
tion. 

Strictly speaking, the quasiparticle approximation is 
sufficient only for the analysis of lowest-order processes. 
In the model under consideration this includes the tree- 
level decay and tree-level scattering processes, which 
are obtained at the order 0(g 2 ) and 0(X 2 ), respec- 
tively (see Fig.EJ). To ensure the consistency of the de- 
scription beyond leading order one should use the so- 
called extended quasiparticle approximation ffdkW$\ . The 
extended quasiparticle approximation for the complex 
scalar field would allow us, for instance, to describe 
tpb — > b — > tpb scattering processes, which are of order 
0(g 4 ). In this paper we are primarily interested in the 
three-loop vertex diagram, whose contribution to the ef- 
fective action contains the fourth power of the coupling 
g. Using the quasiparticle approximation for the propa- 
gators in the vertex diagram thus induces contributions 
which are already of order 0(g 4 ). Consequently, the ex- 
tended quasiparticle approximation would additionally 
induce contributions of higher order in the coupling con- 
stants. Therefore, the calculation of the leading order 
contribution to the vertex CP-violating parameter does 
not require us to go beyond the quasiparticle approxima- 
tion. 

As has been argued in (30l |. in the same approxima- 
tion one can also neglect the Poisson brackets in the 
kinetic equations, which physically corresponds to the 
Stosszahlansatz of Boltzmann. This leads to a simple 
kinetic equation for the spectral function: 



p a V a D p {X, P )=0. 



(B42) 



Let us note, that the quasiparticle approximation for the 
spectral function (|B41[) is consistent with (|B42[) . 

Neglecting the Poisson brackets on the right-hand side 
of (|B37j) . we obtain the Boltzmann equation for the sta- 
tistical propagator: 

p a V a D F (X,p) = ±[D F (X,p)Z p (X,p) 

-E F (X,p)DJX,p)]. (B43) 



Motivated by the fluctuation-dissipation relation (|B40[) . 
we trade the statistical propagator for the one-particle 
number density: 



D F (X,p)= [f(X,p) + \] D p (X,p) 



(B44) 



In view of (|B42[) . we can then rewrite (|B43[) as an equa- 
tion for the phase-space distribution function f(X,p): 

[p a V a f(X lP )}D p (X,p) (B45) 
= i[^ > (X,p)D < (X,p) - D > (X,p)X < (X,p)} , 

where we have introduced 



D>(X,p) 



D F {X,p)±\D p {X,p) 



(B46) 



Equation (|B45|) is very similar to the Boltzmann equation 
for a real scalar field (30l |. There is, however, an impor- 
tant difference. For negative values of po the distribution 
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function / describes antiparticles: 

f(X, -p) = -[f(X,p) + l}. 



(B47) 



In other words, Eq. (|B45[) describes the time evolution 
of both particles and antiparticles. One can obtain an 
explicit equation for / by changing the sign of po : 



[p a V a f(X,p)}D p (X,p) 



(B48) 



= l[t > (X lP )D < (X,p) - D^X^t^X^)] , 

where we have introduced Ti^(X,p) = Eg(X, — p) and 
taken into account that in the quasiparticle approxima- 
tion D p (X,-p) = -D p (X,p). 



Appendix C: Calculation of the self-energies 

The 2PI effective action is given by the sum of all 
2PI diagrams with vertices as given by the interaction 
Lagrangian and internal lines representing the complete 
connected propagators [38[. The structure of the terms 
of the effective action can be read off the diagrams in 

Fig.ru 



ill 40 — l\J x D 2 (x,x), 

ilf> =-±\ 2 f xy D 2 (x,y)D 2 (y,x), 

ilf 5 = - 3 9m9*J xy G mn (x, y)D 2 (x, y) 
- 1 I g* m g n J xy G mn (x,y)D 2 (y,x). 



(Cla) 
(Clb) 

(Clc) 



iT 



(d) _1 



lg t g j9 * m g:fGV(x, y )G mn (v,u) 



X D(y,v)D(x,v)D(y,u)D(x,u), (Cld) 
where, to shorten the notation, we have introduced 

The self-energies of the complex scalar field are obtained 
by functional differentiation of the effective action with 
respect to the two-point correlation function: 



ST 2 [D,G] 
SD(y,x) 



(C2) 



Differentiating the individual contributions to the effec- 
tive action, we obtain 

Y,( a \x,y) = -iS 9 (x,y)\D(x,x), (C3a) 

£ (6 W) = -^ 2 D 2 (x,y)D(y,x), (C3b) 

^ c) (x,y) = -g ig *G^(y,x)D(y,x) , (C3c) 
& d \x,y) = gi g 3 g* m g;J vu G mn (x,v)G^(y,u) 

x D(y,v)D(u,v)D{u,x) . (C3d) 

The components of the self-energy of the system of real 
scalar fields arc obtained upon functional differentiation 



of the effective action with respect to the components of 
the correlation function: 



IL ij (x,y) = 2i 



ST2[D,G] 
SG^(y,x) 



(C4) 



The result of the differentiation reads 
rO) 



V) = -\gig*D\x, y) - y* 9j D 2 (y 7 x) , (C5a) 

^^y)-hS vu G mn {v,u) 

x [ gtg 3 9* m g* n D(x, v)D(x, u)D(y, v)D(y, u) 
+9*9j9m9nD(v,x)D(u,x)D(v,y)D(u,y)] . (C5b) 

The next step is to derive the spectral and statistical 
components of the self-energies (|C3j) and (|C5j) . Upon 
use of the decomposition (|B14[) and of the analogous de- 
composition of the propagators of the real scalar field, 
one easily obtains a corresponding decomposition of the 
self-energies E^, £' c ) and IT( C ) into the statistical and 
spectral components. Linear combinations of the result- 
ing expressions are presented in Eqs. (|lla|) . (jllb|) . and 
p2ap . The calculation of the spectral and statistical com- 
ponents of and H^ d \ which contain two integrations 
over space-time, is more involved (see also (39[). Decom- 
posing the two-point correlation functions in Eq. (|C3cl|) 
into the statistical and spectral components, we get 32 
terms. Each of the terms must be integrated over the 
closed time path C, see Fig.[TDJ It is helpful to use rela- 
tions like the following, 



du° sign c (x° - u°) sign c (u° - y°) D p (x, u)D p (u, y) 



2sign c (x° -y°) / du° D p (x, u)D p (u, y) . (C6) 



One then finds that ten terms vanish upon integration 
over the contour: one term which does not contain any 
sign c functions; five terms which contain only one sign c 
function; two terms which contain a product of two sign c 
functions both depending only on one of the integration 
variables, u or v; and finally two terms which contain a 
product of three sign c functions but depend only on one 
of the "external" arguments, x or y. For the remaining 
terms integration over the contour C reduces to a "single 
time" integration over a combination of the six regions 
in Fig.[TT] Note that the upper limits of the integration 
never exceed the largest time argument (x° in Fig. ll 1|) 
which ensures the causality of the Kadanoff-Baym equa- 
tions. In most cases integration over a part of the uv 
plane can be easily represented as integration over the 
whole plane, < u, v < oo, if two of the spectral func- 
tions are replaced by the corresponding retarded and (or) 
advanced propagators. There are, however, two excep- 
tions: if the resulting integral is only over region (e) or 
region (/) in Fig.[TTJ Using the identities 



1(e) 

hi) 



/(o+c+e) + 1(b) I(a+b+c) ' (C7a) 
f(b+d+f) + 1(a) ~ I(a+b+d) ' 

(C7b) 
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(c) / 
/ (f) 




(a) / 
/ (b) 


(d) 




y x u 
The integration plane in the case x° > y° 



and the definitions of the retarded and advanced propa- 
gators, we can represent the corresponding contributions 
as combinations of integrals over the whole uv plane. Col- 
lecting all the terms, we obtain expressions presented in 
Eqs. pic)) and (pb]) . 



Appendix D: Wigner transformation 

To calculate the self-energies entering the Boltzmann 
equations, we have to Wigner transform products of sev- 
eral two-point functions. Using the definitions (|B24|) . we 
obtain for the Wigner transform of a product of n func- 
tions of the same arguments [3(3] : 

h (x, y)... f n (x, y)^J dll pi . . . dn Pn (2tt) 4 

xS'b-JH-... Pn)f(X, Pl ) . . . f(X, Pn ) . (Dl) 



Equation (|D1[) allows us to Wigner transform the self- 
energies (fTTajl . (fTTbl) and (jl2"a)l . The self-energy (fTTcl) 
has a more complicated structure: 

f(x,y) = I vu h{y,'")f2(u,v)f 3 (u,x)U(y,u)f 5 (x,v) . 

(D2) 



We will now calculate the Wigner transform of (|D2|) in 
the Boltzmann approximation. That is, in each /„ we 
will neglect the deviation of the corresponding center co- 
ordinate from X = X xy . For instance: 

)^fi(X xy ,s yv ). (D3) 

In this approximation the integration over u and v in- 
duces two conditions on the momenta: p v = p u = 0, 
where p u — p 2 + P3 - Pi and p„ = pi + p 2 + P5- In - 
tegration over the relative coordinate s, see Eq. ()B24a[) . 
induces an additional constraint: p = p Sl where p s = 
\{P5 ~ Pi — P3 — Pi)- Thus, in the Boltzmann approxi- 
mation the Wigner-transform of (|D2[) takes the form: 



f(X,p) = J dU pi . . . dn p5 (2^) 4 ^(p„)(2^) 4 ^(p„) 

x (2n)H^p- Ps )f 1 {X,p 1 )... f 5 (X, P5 ) (D4) 



As far as decays are concerned, two of the momenta in 
(|D4|) correspond to the initial and final states, whereas 
three of the momenta correspond to the internal lines of 
the loop. The Dirac 8 functions in (|D4j) ensure conserva- 
tion of four-momentum in each vertex of the loop. 
The self-energy (|12b[) has the structure 



f vu fi{v,u)f 2 (x,v)f 3 (x,u)U(y,v)f 5 (y,u) . 

(D5) 



f(x,y) 



Proceeding in the same way, we again obtain (|D4[) but 
now with p v = pi - p 2 - Pi, p u — Pi + P3 + P5 and 
p s = \{p 2 +P3 ^ Pi — P:-,)- This completes the calculation 
of the Wigner transforms of the self-energies. 



Appendix E: Kinematics of the decay 

In this appendix, we discuss the kinematics of the de- 
cay tpi — > bb in the rest (in the early universe - comoving) 
frame of the thermal bath and in the rest frame of the 
decaying heavy particle. For simplicity, we assume that 
the masses of the toy baryons are negligibly small in com- 
parison with the mass of the heavy real scalar. Let us 
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FIG. 12: Decay of the heavy scalar in the rest frame of the 
plasma and in the rest frame of the decaying particle (see also 
Fig. H] for the assignment of the different momenta). 

first consider the decay in the rest frame of the thermal 
bath. Energy-momentum conservation tells us that the 
momentum p\ of the decaying particle and the momenta 
P2, P3 of the decay products must lie in the same plane. 
The latter depend on the scattering angle see Fig.[T2k. 
Denoting the angle between the bisectrix of the scatter- 
ing angle and the momentum pi of the decaying particle 
by 6 we obtain 



| i^— cos o 



(El) 



The angles between the momentum of the decaying par- 
ticle and those of the decay products are given by a = 
i)/2 + 8 and p = <&/2 - 5. If S = then a = p. If 
6 — > ± \ , then a — > tt and P — > and vice versa. Energy- 
momentum conservation implies, that the energies of the 
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decay product are related to the angle 5 by 



-^2,3 = 



Ex 



IT 



|pi| sin (5 



v/^-|pi| 2 cos 2 5 



(E2) 



If pi = 0, then the energy is equally distributed between 
the decay products. This is also the case if pi ^ but S — 
0. In any other case the energy is distributed unequally. 
In particular if 8 ?, then E 2y s = h(Ei =p |pi|), so that 
in the ultrarelativistic limit one of the decay products has 
almost zero energy, whereas the other receives almost all 
energy of the decaying particle. 

As follows from Eq. (|2l>|) . to calculate the CP-violating 
parameter, we need to evaluate the distribution functions 
associated with the statistical propagators contributing 
to the vertex loop. As discussed in Sca MI) only those two 
terms of (f2"5|) for which both intermediate toy baryons are 
on-shell, whereas the real scalar is off-shell, contribute 
to the CP- violating parameter (see Fig.Hk). Since both 
momenta k\ and fc 2 of the internal toy baryon lines are 
on-shell, the kinematics of these intermediate states is 
the same as the kinematics of the decay products. How- 
ever, due to the presence of the intermediate off-shell real 
scalar, the corresponding scattering angle does not need 
to be equal to fl and the angle between the two scattering 
planes can differ from zero. It is somewhat easier to per- 
form the simultaneous analysis of the intermediate and 
final states' kinematics in the rest frame of the decaying 
particle. The Lorcntz transformation between the two 
frames reads 



A 



J_( E 1 -|Pi| 
M \ -|pi| E x 



(E3) 



As follows from Eq. (|E1[) . in the new frame the scatter- 
ing angle is = tt for both intermediate and final states, 
whereas Eq. (|E2j) implies that the energies are equal to 
M/2. Using the fact that components of the momentum 
orthogonal to the direction of the boost are invariant un- 
der transformation (|E3[) . we can calculate the angle 8 in 
the new frame: 



;| P3 |-|p 2 |)/|pi|. (E4) 



Ei sin 6 

smd = — ; 

y/Ef - |pi| 2 cos 2 <5 



It then follows that 5' > 6 since the denominator of (|E4[) 
is smaller than Ei. 

For a homogeneous and isotropic system the one- 
particle distribution functions depend only on the 
Lorentz-invariant product ku, where k is the particles' 
momentum and u is the four-velocity of the thermal 
bath's rest frame with respect to the chosen frame of 
reference. In particular in thermal equilibrium: 



r(k) = [exp((ku-fx)/T)-l]- 



(E5) 



In the rest frame of the gas u = (1,0,0,0), and we re- 
cover the usual Bosc-Einstcin distribution. Applying the 



Lorentz transformations (|E3[) . we can deduce u in the rest 
frame of the decaying particle 



u = M- 1 (E 1 ,-pi), 



(E6) 



Introducing an orthogonal coordinate system, as is de- 
picted in Fig. 1121 we can then express the arguments of 
the distribution functions in the form: 

uki^ = \ [Ei + |pi | (sin 9 cos <p cos 5' =p cos 6 sin 6')], (E7) 

where tp is the angle between the scattering planes of the 
intermediate and final states (not depicted in Fig. [T2"|) . 



Appendix F: Thermal average of the CP-violating 
parameter 

The CP-violating parameter in vacuum is momentum- 
independent due to Lorentz invariance. Since the sur- 
rounding medium defines a preferred frame (its rest 
frame), the effective CP- violating parameter in medium 
depends explicitly on the momenta of the participating 
particles (see Fig. 2]). In order to investigate the order 
of magnitude of the medium corrections, we consider the 
thermally averaged CP-violating parameter: 



/ effl^ dlil 2 dH| 3 w (pi ,p2,p 3 )e i {p 1 ,p 2 ) 



(Fl) 



/ dU^ dn% 2 dH| 3 w(pi , p 2 , p 3 ) 

where e%(pi,P2) = ^ ac + SeY ied (pi,P2) and w represents 
the gain or loss term (they are equal in equilibrium) . For 
the decay processes: 

w{p 1 ,p 2 ,P3) = {2n) i 6(Ei - E 2 - E 3 )S(p 1 - p 2 - p 3 ) 

x/( P2 )/(P3)[l + /^(pi)]. (F2) 

In the hierarchical limit Mi <IC M2, the CP-violating 
parameter ei only depends on |pi|, see Eq. (|31[). Fur- 
thermore, we set m — > here. Integrating in (|F1 1) over 
momenta of the final states, we find in this approxima- 
tion: 



J °° dq u{q)ei{q) 
J °° dq uj{q) 



where q = \pi\- To calculate the weighting function oj(q) 
we insert thermal distributions characterized by a com- 
mon temperature T and zero chemical potential. We 
consider two cases: 

(i) in the first case we insert Bosc-Einstcin (BE) dis- 
tributions for / and in (|F2|) ; 

(ii) in the second case we insert Maxwell-Boltzmann 
(MB) distribution for / and neglect fy. in (|F2[) . 

The resulting expressions for the weighting function read: 

W (q)= { E t nUh>{®)^ \ si nh(^)J ' ( F3 ) 

MB , 



q 2 exp 



T 



where E q = \/M\ 
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c). The total number density is approximately conserved. 



Appendix G: Numerical details 



To solve the system of Boltzmann equations we 
introduce the transformed variables x = a(t) and ki — 
Sa(t)\pi\ for time and momentum. The constant factor 
S is chosen such that Sx = T^ 1 = (2 • 1.66^/M P it)5. 
The distributions as functions of the transformed mo- 
menta are then well represented, in some sense, in 
the range ki ~ 0.025 — 50.0. 17 In addition, we in- 
troduce the transformed on-shell energies and masses, 
to,; = SxM, and fc? = (k 2 + m, 2 )^ = 5x(|pi| 2 + 
M 2 ) 2. In these coordinates the Liouville operator for 
Robertson- Walker space-time takes the form L[/](|p|) — ?> 
S^ 1 Hk ( l df(ki)/dx, where f(k\) is the transformed one- 
particle distribution function dependent on k\ and x. 
Defining L[f](ki) = df(ki)/dx, the Boltzmann equations 
can be written in the form L[f](ki) = C[.f.](ki) with 
transformed collision integral C[.f.\. 



Homogeneity and isotropy can be exploited to simplify 
the collision integrals significantly. In [7(|, it has been 
shown how the various collision terms for decays, inverse 
decays and 2 — 2 scattering can be reduced to lower di- 
mensional integrals in general. Here, we transform the 
integrals to the new coordinates at the same time. In 
particular, the collision integrals for a scattering process 
12 O 34 (here bb O bb, bb O bb and bb O bb) can be 



17 In particular we require that the approximate numerical value 
of the moments l|G7j l are close to their true values for close-to- 
cquilibrium distributions. Also we demand that particles created 
in decays are not produced with momenta outside of this range 
to a significant extent so that total number densities show the 
expected behavior. 
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FIG. 15: Energy densities of the various species as functions 
of Mi/T for n ~ 0.366 (case c). The ratio of the total energy 
density pt + pi + p^ and the total cosmological energy density 
p is not constant. This feature is due to the different scaling 
behavior of relativistic and nonrelativistic species. For this 
reason the ratio p^,/p increases slightly before the particles 
start to decay. This is more pronounced for smaller washout 
factors (see Fig. 1191 1, 



reduced to a twofold integral: 

CW>34 [./l.Pl) = 

1 1 f f k 3 dk 3 kjdkj n(j Q \ n v 

SHv* 64ir 3 k<{ J J "fe^ k*~ U \ k, 2 ~ m 2j^l2^34X 

{[1 + A(fci)][l + h{k2)]h{h)~h{ki)- 

A(fci)/ 2 (fc 2 )[i + / 3 (fc 3 )][i + h(h)}} , (gi) 

where fc§ = fc§ + k% - fcj and k 2 = [{kf) 2 - m 2 2 ]^- . 
The integrated scattering kernel for a constant 

(momentum-independent) amplitude A and for masslcss 
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FIG. 16: Entropy densities of the various species and the 
total entropy density (s;, + sg + s^,)/s as functions of Mi/T 
for k ~ 0.366 (case c). 
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FIG. 17: Deviation of the distribution function from equi- 
librium for washout factor k ~ 0.366 (case c). 



species 1, 2, 3 and 4 is given by 

D 12 ^z i (k 1 ,k 2 ,k ? , 1 k i ) = 2^0(k 3 + k 4 - \ki -k 2 \)x 

0(h + k 2 - |*3 - fc 4 | )(fc 3 + fc 4 - |fc 3 - feil - |fc 4 - fci|) . 

(G2) 

Similarly, the collision integrals for a particle created in 
inverse decays, 1 f> 23 (here ipi -H- bb and ifti -O- 66), can 
be reduced to a single integral: 

CW23[./i.](*i)=# saSEf/^^-majDi^asX 

{[l + A(fci)]/2(fc 2 )/3(M- 

/i(fci)[l + /2(fc 2 )][l + / 3 (fc 3 )]}, (G3) 



FIG. 18: Number densities of the various species and the 
generated asymmetry rj as functions of Mi/T for k ~ 0.01 
(case a). 
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FIG. 19: Energy densities of the various species and the total 
energy density (p b + Pi + Pi>) / P as functions of Mi/T for 
k ~ 0.01 (case a). 



where fc° = ^i — ^3 an( l ^2 = \{k 2 ) 2 — m 2 2 ] 2 . The inte- 
grated scattering kernel Di^ 23 is given by 



Di^ 23 (ki,k 2 , k 3 ) 
2A 



i-\k2-k 3 \)e((k 2 + k 3 )-h) 



(G4) 



Finally, the collision integrals for a particle created in 
decays, 12 f> 3 (here bb -o- -0i and bb O "01 ), can be 
reduced to the single integral 

Ci2« 3 [./i.](fci) =#3^0 / ^f 1 ^ - m 2 )£> 12 ^ 3 x 

{[l + /i(Ai)][l + / 2 (*a)]/3(A3)- 
/i(fci)/2(k)[l + /3(*s)]}, (G5) 

where fc° = ^3 — ^1 an( l ^2 = \{k 2 ) 2 — m 2 2 ]5. The inte- 
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Sec. lIVl we assume that the interactions 66 O 66, 66 «-» 66 
and 66 •H- 66 are rapid enough to keep the distribution 
functions of 6 and 6 very close to their equilibrium dis- 
tributions, parametrized by ao,ai,ao and a\: 



/ a e ?(fci) = [exp(a 
/°*(fa) = [exp(d 



■ aiki) — 1] 1 
•Ox&i) - 



(G8) 



Assuming that 66 -O- 66 alone is much faster than the in- 
verse decays into tpi , the evolution of / and / can there- 
fore be described by means of three parameters ag , do and 
a±. The equations for the evolution of these parameters 
are obtained by forming the moments n[.] of Eqs. (|36aj) 
and (l36bl): 18 



FIG. 20: Entropy densities of the various species and the 
total entropy density (s& + s b + s^)/s as functions of M\/T 
for k ~ 0.01 (case a). 
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FIG. 21: Deviation of the distribution function f^ >1 from equi- 
librium for washout factor « ~ 0.01 (case a). 



grated scattering kernel Z?i2-<->3 is given by 



A.2-B-3(fa,fa, fa) 

2.4 



fa 



(fa-|fa-fa|)e((fa + fa)-fa) 



(G6) 



Number density and energy density corresponding to the 
distribution / in transformed coordinates read 



«[/] = 2M^) 3 /(fa) 2 />i)*i, 
p[/H2^(^) 4 /(fa) 2 fci/(fa)*i 



(G7) 



For massless particles these are the second and third mo- 
ment of the distribution, respectively. As outlined in 



dx 


. 9a _ 


dai 

+ —^n 
dx 


dai 


n 










dao 
dx 


~dfl q 


dai 
dx 


~dfl q 
dai 


n 











Here, we used 



a 



bb^bb 



[fJ] 



(G9) 



Cbb^bbif] = n C bb ^ bb [f] = and 



'bb++bb[fJ}\ = °- The third equa- 
tion is obtained by forming the moment /?[.] of the sum 
of Eqs. QUID and (p6b| . i.e. 



P [Ufa 11 } 


+ P 


m q }_ = 


dao 
dx 


\dfa q ] 

da . 


da Q 
dx 


dai 
dx 


~dfl q 
_ dai 


dai 
dx 


dai 







eg 



dao 



(G10) 



where we used p C bb ^ bb [f , f] + p C bb ^ bb [f,f] = 0. 

The derivatives of f^ q with respect to the parameters 
can be rewritten as 

df e a q (ki) 



da,; 



(fci ) r[l + /r(fc 1 )]/r(fc 1 ),i = 0,l. (Gil) 



An analogous relation holds for the derivatives of f^ q 
with respect to do and a\. Solving Eqs. (|G9|I and ()G10p 
for dao/dx, ddo/dx and dai/dx, we find the differential 



18 Here and in the following we use the abbreviations Cf >b ^.^ 1 = 
C bb ^ i , 1 [fa q ,/VJ and Css^^j = C u ^ l jfyj. Also note 
that /, / and /y, are functions of the transformed coordinates, 



here. 
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equations for the three parameters: 

dao ai4fc?(l + iT)jT] + n[CW v , 1 ] 
^ ~ n[(l + f?)f«\ 

dao ain[fc°(l + /^)/ a ^] +n[C^. J 
dx ' n[(l + 

^ = -((„[<% 5 ^ 1 ]p / + n[A?(l + / a e9 )/ a e9 ]pc) 

+ (4^v>>/ + n[*?(l + ^)/^]pc) 

x n[(l + / a e ")/ a e "]p[(l + /l 9 )/^]) A + Po/P/, (G12) 

where we have defined 

h = Pf (n[(i + r a q )r a q H(i + f?)f?]p f 

+ n[k° 1 (l + f?)f?]n[(l + f a i)r a i] P [(l + %>)%>] 

+ n[fc?(i + .o/ a e9 ]n[(i + r a q )r a q ]p[(\ + r a q )r a q ]) , 

(G13) 

as well as 

Pf = + f?)f?] - p[kl{l + j?)/?] , 

Pc =p[C bb +ni, 1 ] + p[C- b i^ itl ] . (G14) 

As stated in the main text, we need to start with finite 
chemical potentials as to avoid the occurrence of Bosc- 
Einstcin condensation. We choose the minimal accept- 
able value a = ao = 0.5, corresponding to p b = p b = 
—0.5 To and p^ 1 = 2p b . The initial value a\ = 1 corre- 
sponds to the initial cosmological temperature To. We 
checked that the results do not depend on T as long as 
T) ^> M\. The heavy species %pi is subject to relatively 
weak interactions only, so that its distribution function 
can deviate from kinetic equilibrium. Therefore, we solve 
the full Boltzmann equation for ipi, 

HUiKki) = C^t+bblffa, f](ki) + C^^nlffa, f](ki) , 

(G15) 

along with the integrated ones for b and 6. 

Because of the integration of the equations for the 
massless species all collision terms for 2 — 2 scattering 
drop out of the system. In order to verify that the rates 
for these processes are much larger than the ones of the 
decays and inverse decays we have computed the rates 
for these processes numerically. The maximum (during 
the full evolution) of the ratio of T^^^ and T bb ^ bb is 
exemplarily presented in Fig.Q2] (the rates for the other 



2 — 2 processes are similar). It shows that we can choose 
A ~ 1 or smaller for most of the relevant range of \g\ 2 if 
we demand that r bb -++ bb -/Tbb++ii> 1 ^ 10 3 as criterion that 
b and b arc in kinetic equilibrium at all times. Here the 
equilibrium shape of / and / is not distorted by the ex- 
pansion since we are dealing with massless particles. In 
addition, it can be argued that the 2 — 2 processes are 
meant to model rapid gauge interactions with different 
particles which would have the same effect of equilibrat- 
ing b and b. In this sense wc could even formally tolerate 
nonperturbative values of A. 

To turn the equations into a system of ordinary differ- 
ential equations (ODE) the distribution functions were 
discrctized on a grid of dimension 400 with linearly 
increasing spacings in the range fci ~ 0.025 . . . 50.0 
to account for the characteristic behavior of close-to- 
equilibrium distributions at small and large momenta. 
All integrals were approximated by Riemann sums on this 
grid. The system of Boltzmann equations behaves nu- 
merically stiff. This means that it is advisable to use an 
implicit method for its numerical solution to achieve ac- 
ceptable step sizes (and hence acceptable execution times 
and numerical errors). Here CVODE with its backward 
differentiation formula with Newton iteration was used 
as ODE solver. The full Jacobian was computed ana- 
lytically in every external step. A relative tolerance of 
10~ 8 was attributed to every momentum mode. Because 
of the implicit method all solutions were computed in 
£>(10 3 ) steps. 

Since the global systematic error due to the discretiza- 
tion cannot be computed within the method the proper 
behavior of the system was tested by successive refine- 
ment of the grid and comparison of some of the macro- 
scopic quantities with the theory predictions. For this 
purpose, we present two examples of the number densi- 
ties n x , the energy densities p x and the entropy densities 
for the washout factors k ~ 0.366 (case c in Fig. fT4lfr6^ 
and k ~ 0.01 (case a in Fig. [T8ll2"0"|) . The total number 
density (n b + 2n 1 p l + n b )/s is almost conserved (as dis- 
cussed in Sec. lIVp . The ratio (p b + p^ 11 + p b )/p is not 
constant (see Fig. ITS"!) . This behavior is expected for a 
system involving nonrelativistic massive particles and is 
also observed for the bottom-up equations. The ratio is 
much smaller than one so that it is justified to neglect 
the backreaction on the curvature. Finally, the total en- 
tropy density is steadily increasing as it should. Figures 
[T7l and [2T1 show the deviation of the distribution function 
/ 01 from kinetic equilibrium ones for which the curves 
would be straight lines. The deviation from equilibrium 
is larger for smaller values of n and increases at late times, 
as expected. 
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